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Chapter 1 

Introduction. 



In the last decades, new ideas and tools coming from quantum informa- 
tion community [BDOO, NCOO, :2001, ERD03] have brought new insights and 
a complementary point of view to problems in quantum physics, in general, 
and highly correlated quantum systems [OAFF02, ON02], in particular. In 
fact, it was point by Preskill[PreOO] that this interdisciplinary area of re- 
search could bring better understanding to questions that appear in many 
body quantum entanglement and deeper classifications of different phases 
that appear in strongly entangled systems [Lau83, Wen02]. 

Following these lines, the initial purpose of this work was to study entan- 
glement properties in a many-body quantum system characterized by some 
local interactions. For this aim, we studied low dimensional magnetic models 
that can be described by quantum spin chains. Their dynamics are defined 
by some parameters that, depending on their values, can drastically change 
the behavior of the ground state and their correlations. In fact, these mod- 
els, in spite of their simplicity, display a rich structure in the phase diagram, 
e.g. ordered-disordered magnetic phases or phase transition due to the for- 
mation of vortices. Correlations become long range at the critical points and 
when the temperature of the system is fixed to zero, only pure quantum cor- 
relation, i.e. entanglement, can drive the system to the phase transition. So, 
entanglement is a crucial characteristic in the behavior of quantum systems. 
We show that the phase diagram of a given system can be described just 
using their entanglement properties. Due to the complete description of the 
correlation in one dimensional model, we could also explained the success 
of some numerical variational method, as density matrix renormalization 
group, in simulating these systems. The core of this chapter is contained in 
the published works: "Entanglement in quantum critical phenomena" by G. 
Vidal, J. I. Latorre, E. Rico, A. Kitaev[VLRK03], "Ground state entangle- 
ment in quantum spin chains" by J. I. Latorre, E. Rico, G. Vidal[LRV04] 
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and "Entanlement entropy in the Lipkin-Meshkov-Glick model" by J. I. La- 
ter-re, R. Orus, E. Rico, J. Vidal[LORV04]. 

In addition to the first point, special mathematical tools, developed in 
the area of field theory, can be applied to the study of critical phenomena 
in one dimensional quantum systems. In our works, we link conformal field 
theory concepts with entanglement behavior. Conformal field theories have 
succeeded to describe any property of almost any critical model from the 
knowledge of general symmetries of the system. As it will be shown, entan- 
glement does not escape from this classification and so, it is also related with 
the deep ideas of symmetry and universality. Another tool that it is espe- 
cially important in the description of phase transition is the renormalization 
group. We proposed how renormalization group flows can be characterized 
with tools of quantum information theory as majorization and entropy. In 
fact, it seems that using concepts like majorization, it can be achieved a 
better understanding in the flows of the hamiltonians or how low energy 
properties of a system appears as we change the scale (energy, momenta or 
length) when we prove it. The main part of this chapter is described in: 
"Fine-grained entanglement loss along renormalization group flows" by J.I. 
Latorre, C.A. Lutken, E. Rico, G. Vidal[LLRV05]. 

The last chapter is focused in the analysis and study of matrix product 
states which is a formalism to represent the state of any quantum system. 
This method is specially well suited to describe in an efficient way one di- 
mensional quantum systems with translational symmetry. In particular, we 
will address the possibility to realize renormalization group transformation 
at the quantum level on the state that describes the system, independently 
of any dynamics. Also, we will show how the fixed point of the transfor- 
mation can be characterized and we will give a complete classification of 
quantum states that are fixed point in some relevant cases. The most im- 
portant aspects of this chapter appears in the article: "Renormalization 
group transformations on quantum states" by F. Verstraete, J.I. Cirac, J.I. 
Latorre, E. Rico, M.M. Wolf[VCL+05]. 



Chapter 2 



Scaling of entanglement in 
(l+l)-dimensional systems. 



2.1 XXZ and XY models. 

The analysis of entanglement properties in the ground state of two ex- 
actly solvable models, XY and XXZ models, is one of the main aims of 
this work[VLRK03, LRV04]. These models describe quantum spin chains 
with first neighbor interactions and an applied static magnetic field in the 
z-direction. Their hamiltonians can be written as, 

N 



j=l ^ / 

N (\ A 
Hxxz = E U + + 7°?°?+i " 2 ° 

3=1 V 



(2.1) 



where {(J ( j} a =x,y,z are the Pauli matrices, such that [ct°,ct^] = 2i5jke a ^ e aj. 

This implies commutation relations in different sites [aj, ] = with j k 

but they anticommute at the same site {a",^} = 2I5 a @ . The parameter 
7 quantifies the anisotropy in the interaction and A the strength of the 
magnetic field, and we are considering periodic boundary conditions, i.e. 

2.1.1 Bethe-Ansatz in the XXZ model. 

There are several techniques to determine the low energy behavior of the 
XXZ model. Recent mathematical tools, like bosonization and mapping 
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2.1. XXZ and XY models. 



to sigma models (see for example [Pra91, Aue94, GNT98] and references 
therein), have shown the continuum theory that describes this model. Nev- 
ertheless, Bethe[Bet31], in 1931, was the first to show how to get the ther- 
modynamical properties using the symmetries that characterize the model 
and that make it integrable. The hamiltonian of the model can be split as 



where Hxx corresponds to the interaction in the xy-plane of the spin space 
and it is the limit of zero anisotropy and magnetic field, i.e. {7, A} — > {0, 0}; 
Hising describes the interaction in the z-direction of the spin and it is the 
main contribution in a highly anisotropic model, i.e. 7 > 1; finally M z 
corresponds to the total magnetization in the z-axis. 

In the XXZ model, the hamiltonian, H(k) = Hxx + K H K , varies as a 
function of a dimensionless parameter k (the magnetic field or the anisotropy) ,| 
where Hxx and H K commute, so they can be simultaneously diagonalised. 
While the coupling k is modified, the energy of states changes. When a 
first excited level decreases its energy enough to get the ground state value, 
the two levels cross each other (fig. 2.1) and a point of non-analyticity in 
the ground state is created as a function of k. After and before the level- 
crossing, the ground state is an eigenstate of some properties like the mag- 
netization but this is no longer true at this specific situation. Any point of 
non-analyticity in the system is identified as a quantum phase transition. 

Figure 2.1: Sequences of level crossing in the XXX model, 7 = 1, while 
the magnetic field increases. The initial points in the vertical axes indicates 
the different spin sectors in which the spectrum of the theory can be orga- 
nized. The different lines that appear from them indicate the splitting of 
the degeneracy due to the Zeeman term. 



As we have just seen, the total spin angular momentum in the z-direction, 
M z = 7} 5Zj=i ° Z ji 1S a good quantum number to block-diagonalize the inter- 
action. The second symmetry that is explicitly used in the Bethe Ansatz is 
the translational invariance of the system by any number of lattice spacing 
and the periodic boundary conditions. 

Due to these symmetries, any translational invariant eigenstate of the 
total magnetization in the z-direction can be an eigenstate of the XXZ 
model, so that, the ferromagnetic state, i.e. \F) = \ ]] ... |), should be 
an eigenstate of the hamiltonian. In fact, it can be shown that H\F) = 
(^f ~ A f ) \F) = E F \F), so it is an energy eigenstate, and M Z \F) = f \F), 




(2.2) 



Hxx + lH Ising - \M Z 
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it has maximum magnetization. From the ferromagnetic state, it can be 
derived a translational invariant state with a well defined magnetization 
applying the operator a~ that decreases the magnetization in one unit, 

1 N 1 N 27TT1 

L£ e *V|F> = 4=£^ * = = ne{l,N} (2.3) 



1=1 v i=i 

where \l) is the state with the spin at the position / reversed respect to the 
ferromagnetic state and with energy (E — Ep) = — (7 — cos k) + A. In the 
case with two reversed spins, the wave function appears as, 

|tt>= <hMh,h)= £ a{hM)^T 2 \ F ) ( 2 - 4 ) 

l<h<l 2 <N l<h<h<N 

Bethe's idea was to consider a superposition of plane wave functions: a(li, I2) =| 
Ae^ klh+k21 ^ + A' e^ klh+k2ll \ Imposing that |*) is an eigenvector of the 
hamiltonian, i.e. H\^f) = E\fy), then 



E + 1 










£ + 7 





a(h,h) = a(h,l 2 + l)+a(h - 1,Z 2 ), h = h + 1 

a(h,l 2 ) = a(h,l 2 + 1) + a(h,l 2 ~ l)+ ^ 
+ a(h + l,h) + a(h-l,h), l 2 >h + l 



where E — Ep = 2A — X^j=i 2 (7 ~~ cos ^j)- Applying the Bethe's ansatz into 
the last equations, an equivalent relation is derived 

2ja(h,h + 1) = a{h,h) + a(h + l,h + 1) (2.6) 



or 



A' e i{k\+k2) _|_ 1 _ 2je ik2 6 ' 

In a general case, where there are r reversed spin respect to the ferromagnetic 
state, the eigenvector is written as 

|*)= a(h,...,l r )\h,...,l r ), (2.8) 

l<h<...<l r <N 

where the ansatz for the coefficient in the wave function is 



a(h,...,l r ) = Y exp 1^^ + -^^!, (2.9) 
veSr \ j=i i<j 

V G S r denotes one of the r! permutations of {l,...,r} and ki and Oij 
with € {!,... ,r} are the parameters to be determined. Three general 
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2.1. XXZ and XY models. 



conditions hold for these parameters [Orb58]: 
9i j 7 sin fc ' 9 fcj . 

2 cos ^i + 7COS ^_^ (2.10) 
Nki = 2TTli + ^2e itj i€{l,...,r}, 

where the integers Tj are called Bethe quantum numbers and completely 
determine the state. It is known [YY66] that the set {li} with the lowest 
energy for each M z = N/2 — r satisfies: 

N 

l i = M z -l + 2i = — - r -l + 2t i£{l,..,r}, (2.11) 
and the energy reads 

E = 2_,(7 - cos fci) - \M Z . (2.12) 

i=l 

Therefore, the values 0's and k's give you the whole description of the 
ground state of the XXZ model. In this work, we solved the previous non- 
linear equations system, eq.(2.10), with a genetic algorithm, and the maxi- 
mum numerical error allowed in the parameters 0's and fe's was 1(T 5 . Also, 
the spin chains that we considered had an even number of sites N, so the 
eigenvalue of the M z operator in the ferromagnetic system is the integer 
N/2. 



2.1.2 Diagonalization and continuum theory in the XY model.| 

In this section, we will fix our attention in the XY model, 

The model without magnetic field was diagonalized in 1961 by Lieb, Schultz 
and Mattis[LSM61]. In 1962, Katsura[Kat62] solved it with an applied mag- 
netic field. Pfeuty[Pfe70] solved the Ising limit with a transverse magnetic 
field in 1970. In 1971, Barouch and McCoy[BM71] got the spin correlation 
functions. See also the textbooks [CDS, CH]. 

To diagonalise this hamiltonian several unitary transformations have to 
be applied. Due to the properties of spin operators that mix commuta- 
tion and anticommutation relations, it is particularly useful to describe one 
dimensional spin system by one dimensional spinless fermions, i.e. with 
anticommuting operators. This step is performed with the Jordan- Wigner 
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transformation that maps the Hilbert space of a spin system into a fermionic 



one, 



x ■ y 

II < ) ° l ~ 2 ; {^m,aj} = S ml , {a ; ,a m } = V{Z,?n}. (2.14) 

\m<l / 



Then the hamiltonian is recast into 

x ^ / o +1 Q|+O Q| + i Qj+iQi+Q Q +i A+A \ f . 

H xy = }_^\ ^ + 7 2 Aajaj I . (2.15) 

The Jordan- Wigner transformation changes the boundary conditions de- 
pending on the total z-component of the spin, 

a*+i=( II <)? 1 ^-=i- N e i ^=^a 1 . (2.16) 

\rn<N+l J 

Then, in the sector with ~ Yl m =i °m = 0> an d y even (odd), the system has 
periodic (antiperiodic) boundary conditions. 

Due to translational invariance, applying the Fourier transformation al- 
most diagonalizes the hamiltonian, 

N 

4 = 5>ze"^ {d k ,S p } = 5 kp , {d k J p } = V{k, P }, (2.17) 
i=i 

then, the hamiltonian is recast into 

N/2 



h xy = ( cos ^jf - A ) 44+ 



k=-N/2+l 



N/2 

h sr^ ■ Ink / ~ ~ * * 
} J sin — - [dkd-k + d k d'_ 



(2.18) 



2 

fc=-JV/2+l 



If 7 = 0, the theory describes what is called XX model. This model has 
U(l) invariance which is a continuous symmetry. Coleman-Mermin-Wagner 
theorem shows that a continuous symmetry can not be spontaneously broken 
in one dimensional system. Then, there is no transition from a disordered to 
an ordered phase. Nevertheless, it is still possible to have polynomial decay 
in the correlation functions via a mechanism described by Berezinskii[Ber71], 
Kosterlitz and Thouless[KT73] with the formation of vortices. In this limit 
the hamiltonian is completely diagonal, 



N/2 . . 

Hxx= ( cos ^- A )44, (2.19) 

k=-N/2+l ^ 
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The thermodynamical limit is obtained letting — > (j) and Y^k=-N/2+i \ 

rir dip 
J-n 2tt' 

Hxx = f ^ (cos <j> - A) 4^ = jT ^ (2.20) 

Once the hamiltonian is diagonalized, it is straightforward to obtain the 
continuum limit. The spectrum of the XX model has two Fermi points 
<f>F = iarccosA, i.e. where the energy goes to zero. To develop the 

I 



low energy theory, we define two fermionic modes 
I 



= Rj, and 



= L<j), i.e. right and left movers, and expanding the spectrum 

<£->■- \4f\ 

around the Fermi points, the low energy theory turns out to be 
Hxx — ► 

VF J dct)ct> (fi^fy - = -iv F J dx (h ] d x R - Pd x i) 



with vf = and x the position in the coordinate space. From this 

<t>F 

expression, the XX model can be seen as a free theory of two noninteracting 
fermionic fields and it is massless in the interval |A| < 1. For |A| > 1 the 
spins in ground state are polarized in the direction of the magnetic field. 

For the general case of parameters (A, 7), a Bogoliubov transformation 
recast the theory into a free fermionic one, 

h = cos — d k - 1 sin yol_ fc > 

{bp, Sj.} = 5 pk , {b p , b k ] = V{k,p}, ^ ^ 

cos^-A 

cos o k - 



7(cos^-A) 2 + 7 



2 „; T1 2 2wk_ 



then, up to an overall constant, the hamiltonian for the XY model can be 
rewritten as 



H XY = Yl ^kblk, A fe = l/( cos ^" A ) +7 2 sin 2 ^. (2.23) 



k=-N/2+l 

and in the thermodynamical limit 



H XY = J W ^ A = ^(cos0-A) 2 + 7 2 sin 2 </>. (2.24) 
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Figure 2.2: Energy spectrum for different values of the anisotropy 7 
and magnetic field A in the XY model. The ones with zero energy gap 
correspond to the critical points. In the figure, critical XX model (red line), 
critical Ising model (dashed green line), the case [7 = 0, A = 1] (dot-dashed 
blue line), the Ising case without magnetic field (straight blue line) and the 
case [7 = 0.5, A = 0.5] (dashed violet line). 



For the continuum limit, we consider first the case A = 1, i.e. critical 
XY model. In this case, there is just one Fermi point (pp = 0. Expanding 
the spectrum around this point, the hamiltonian can be written as follows 

Hxy^I J d<\> \4>\ bfa = 7 J d(t> <\> (P^ + bl^b-t) . (2.25) 

Defining two real fermionic modes, i.e. majorana fermions, such that 

J(p>0 V27T V 

a L , x = [ -% (e-^U + J*tL) , (2-26) 
J(p>0 v2vr v r / 

{d a ,x, &f3,y} = $af3 8(x — y), O a ,a;«Q,a: = « • 

the system explicitly appears as a massless majorana theory where the 
anisotropy just redefines its velocity, 

— 27 r 

Hxy -> -y / dx (a R>x d x d RtX - a L:X d x a L:X ) . (2.27) 

For the rest of the parameter space the continuum theory is still recast 
into a relativistic fermionic theory but with a mass term. Expanding the 
spectrum A^, around the point of minimum energy, that is Fermi points, 
it is straightforward to see that two different regions have to be described 
depending on the number of Fermi points, 



JO A + 7 2 > 1 

y ± arccos j^-j A + Y < 1 

For the case 4>p = the low energy theory turns out to be 



(2.28) 



# ^m 2 + 4^ 2 b\b h (2.29) 



where m = 1 — A and uf = A + 7 2 — 1. The system is a massive theory 
with a Klein-Gordon dispersion relation. The mass term describes how 
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much energy has to be given to the system to create an excitation. The 
critical points are those where the mass goes to zero, in this case A — > 1. 
In the case <Pf = ± arccos ^ 2 , two fermionic fields are defined, left and 



right movers, around the two different Fermi points, 



d,i 



= Ra) and 

<M<M 



d 



= La.. Then the low energy hamiltonian is recast into, 



|0F|=arccos — 



(2.30) 



with m = 7 2 ^1 — j^yi^J an d V F = 1 — 7 2 — • ^ n this case, the the- 
ory appears as a massive relativistic model of two noninteracting fermionic 
fields and the critical theory appears when 7^0 where the XX model is 
recovered. 

Figure 2.3: Critical region for the XY model in the plane of the parameters, 
anisotropy 7 and magnetic field A. 



2.2 Ground state and reduced density matrix in a 
free theory. 

When the hamiltonian of a quantum system is diagonalized, it is quite simple 
to obtain the ground state of the system. From eq. (2.23), the XY model can 
be seen as the sum of decoupled fermionic modes, with a positive dispersion 
relation A/% for every mode. Then, the state that minimizes the energy is 
the one that has no population in any of the modes or it is annihilated by 
the 6fc operator, 



h\g) = o, 



b{h\g) 



0, Vfe e 



N N 

T + 1 'T 



(2.31) 



Due to the fact that the system is made of noninteracting fermionic 
modes, the ground state is the tensor product of the occupation of every 
mode, 



\g)(g\ 






1 -n k 



®fce[-f +i,f 



N N' 



(2.32) 
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In fact, the eigenvalues of the density matrix can be obtained from the cor- 
relation matrix of a pair of majorana operators (see [Pes03, Pes04b, KitOO]) 
defined by, 

_ Oirn ~i~ ®m - _ Qm . r ~ - \ r » « _ 1 



1 

&2n = ~m Yl (f^A + f* k b k 

JV 

«2n-l = -7= \hnkb\ ~ h* k b, 



(2.33) 



where / n fc and /i n fc implement the Fourier and Bogoliubov transformations, 
relating the majorana fermions in the coordinate space with the complex 
fermions in the momentum representation and they fulfill, 



JV JV 
2 2 



^ ] fnkfmk ~ y ] h n kh mk — S nm ; 



fc=-f+l k=-f+ 1 

JV JV 

2 2 



^ fnkKnk — ^2 fnkfmk — 9nm', (2.34) 
fe=-f+l fe=-f+l 

JV 

1 i^fn-rn^ //V COS ^ - A - ij sin 3& 



5ram — N ^ 



c »27r(n-m)fc/jV — - /v - / 0111 — 

iV fc= f| +i | cos ^_ A -z 7 sin^|- 

Therefore, the 2N x 2iV correlation matrix of the majorana operators is 
written as follows, 



r lm = ([ x ^ ) ("2m, «2m-l)) 



02n 
fl2n-l 

JV 



JV 



m J 



JV 

2 / j? ^* \ / h* 



1 / fnk f*k ) [ / mfe — ] = 1 f <5 n m ~W\ 

(2.35) 
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2.2. Ground state and reduced density matrix in a free theory. 



Figure 2.4: Steps followed in the text to diagonalize and obtain the eigen- 
values of the density matrix from the correlation matrix. 



In addition to this, because the XY model can be described as a free 
theory, applying the Wick's theorem, the expectation value of any set of 
operators can be decomposed as the sum of contractions of every pair of 
operators, 

{a\a\d 3 a4) = {a\a\) {a^a A ) - (a[a 3 ) {a\a A ) + (a\a±) (a\bz) , (2.36) 

where a ] n = £ fe e ~ i27rkn / N (b[ cos 6 k /2 - ib_ k sm6 k /2} . This means that the 
state is fully characterized by the expectation values of the one body and 
two body correlators, i.e. it is a gaussian state. Obviously, this property 
is still valid although we trace out every degree of freedom except the ones 
considered in the correlation function, i.e. if we use the density matrix 
P1234 = ti'v-i234|5 , )(5 , |- In particular, we are interesting in finite subsystems 
of L contiguous sites out of the TV" sites of the total chain. Again, according 
to Wick's theorem, this property holds if the reduced density matrix of L 
sites is the exponential of a quadratic hamiltonian[CP01], i.e. 



PL 



«- E 



Z 

n,m=l 



(2.37) 



with Z a normalization factor, A an hermitian matrix and B an antisymmet- 
ric one. Because H is a quadratic form, it can be diagonalized in a similar 
way as in the case of the XY model[LSM61], getting a new set of L fermionic 
modes, such that, 



L/2 

U= ^2 ^kb\bk, PL = tr N - L \g)(g\ = Y[p k , n k = tx\dktyj>k 

k=-L/2+l ®fc 



(2.38) 



>kk' 



The real major ana fermions are written in terms of these new fermionic 
modes as follows, 



" 2n = 77f 12 {jnk% + f* nk b k/ 

fc=-4+i 

(2.39) 

02«-i = —J^ zZj (h nk b{ - h* nk b kj 
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and the 2L x 2L correlation matrix of these major ana fermions are related 
with the eigenvalues of the density matrix as in the eq.(2.35), 



2.3 Analysis of entanglement behavior. 

In the last sections, we have seen how to diagonalize the hamiltonian of two 
quantum spin models, the XXZ and XY model, with two different meth- 
ods, the Bethe Ansatz and unitary transformations. These methods allow 
to obtain the ground state of the system for any value of the parameters 
(anisotropy 7 and magnetic field A), that define the interaction, and there- 
fore to obtain a complete knowledge of the behavior of the theory. As 
we have seen, there are some specific values of the couplings that make 
the theory gapless or in other words, the correlation length diverges. This 
divergence is not driven by any thermal fluctuation since in our analysis 
the temperature is always zero and plays no role in these models. Only 
quantum fluctuations in the ground state can lead the system to the crit- 
ical region. This fact is the main reason to describe the parameter space 
of these theories with measures of entanglement or pure quantum correla- 
tions that have been developes in the last decades in the area of Quantum 
Information. For this purpose, given the pure ground state ^(7, A)), we 
perform a bipartition of the system getting the density matrix of L con- 
tiguous spins out of the whole chain, i.e. P( 7> a,l) = ^N-Llgil, A))(<7(7, A)|. 
Then, we employ the entropy of entanglement [BBPS96] as a measure of the 
quantum correlation between both subsystems. This measure is defined by 
e l(\9(7,X))) = S(p il)X ,L)) = -trp( 7jA ,z,) log 2 p (7|Aji ). Changing the number 
L of spins in the subsystem we look for the scaling properties of the entropy. 

Figure 2.5: Entropy of entanglement in a spin chain for a subsystem of L 
contiguous sites. 



Entanglement entropy is a positive defined quantity and it is zero for 
separable state, those one that can be written as the tensor product of its 
parts. In one dimensional systems of N sites that are translational invariant 





(2.40) 
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the entropy of the first L contiguous sites is the same as the entropy of the 
N — L first sites, i.e. Sl = Sn-l- A third property that will appear in our 
analysis is a consequence of the strong subadditivity of the entropy [LR73]: 
entropy is a concave function respect to the number of sites consider in the 
bipartition, i.e. S L > s l-m+s l +m with M £ [ 0)m i n {iV - L,L}}. 

In this section, we will see that just the entropy can describe the complete 
behavior of the system in the parameter space, what allows us to link the 
scaling properties with deeper concepts like symmetry and universality. 

2.3.1 Entropy of entanglement in the XXZ model 

Although Bethe ansatz gives the eigenvalues of thermodynamical quantities 
like energy or magnetization in an easy way, it is not so simple to obtain the 
complete description of the wave function. So that, we were able to compute 
the ground state for systems up to 18 spins. 

In the fig. 2.6, we see how the finite size of the system saturates the 
maximum value of the entropy. The plot shows the curves of the entropy for 
chains with N £ {8, 10, 12, 14, 16, 18} sites. Although in this case, we are 
limited with finite chains, the initial behavior of the curve, where the finite 
size effects have less influence, obeys a logarithm behavior such that 

S L ~±\og 2 L. (2.41) 



Figure 2.6: Entropy of entanglement for an isotropic Heisenberg chain with 
no magnetic field applied for several values of the total number of spins. 
Triangles (N=18), stars (N=16), diamonds (N=14), triangles (N=12), stars 
(N=10), diamonds (N=8) 



Also, we show how the entropy is modified as the anisotropy or mag- 
netic field change. In fig. 2.7, we analyze the XXX model (or Heisenberg 
model) with an applied magnetic field. It is well known that there are two 
limiting cases in this model: when the field is zero the ground state is a sin- 
glet but while the field increases, states with higher magnetization become 
the ground state. When the magnetic field has a value equal to 2, the fer- 
romagnetic state becomes the ground state which is a separable state with 
maximum magnetization. In the thermodynamical limit, the interval |A| < 2 
is a critical interval [BF64], this fact appears in the finite size studying of the 
entropy with a logarithmic scaling in the interval. 

In the fig. 2.8 is represented the behavior of the entanglement entropy in 
an anisotropic spin chain with N = 18 sites and no magnetic field applied. It 
is known [Aff98] that for values |7| < 1, the anisotropy remains as a marginal 
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Figure 2.7: Scaling of the entropy in the Heisenberg model with an 
applied magnetic field. The plots scale with a logarithmic behavior in 
the interval |A| < 2. The upper curve correspond to a zero magnetic 
field while the lower one is for A = 1.97; the other curves correspond to 
A E {0.24,0.68,1.05,1.35,1.59,1.77,1.89}. In the limit case A -> 2, the 
entropy vanishes. 



operator, so the long range behavior of the system does not change while 
the anisotropy is modified. But for | — ^ | > 1, this parameter turns to be a 
relevant deformation and the energy spectrum of the theory develops a gap. 
Again this behavior is described by the scaling of the entropy, for |7| < 1 
the curves overlap with a logarithmic behavior, but for I7I > 1 the curves 
start to bend down and saturate due to the anisotropy. 

Figure 2.8: Entanglement entropy in an anisotropic spin chain with N = 18 
sites. The upper curves that overlap correspond to systems with system 
with anisotropy I7I < 1, while the three curves that bend down correspond 
to 7 = 1.5, 7 = 2.0 and 7 = 2.5 respectively. 



2.3.2 Entropy of entanglement in the XY model 

In this section we will show the scaling properties of the XY model, for which 
it is possible to get the thermodynamical limit as we saw in the last section. 
Several properties make possible a numerical treatment of this model even 
in the infinite length chain. As we have seen in the last section this model is 
described by a free fermion theory and the ground state is a gaussian state, 
only depends on one and two body correlation functions, due to this fact 
the reduced density matrix of L sites is written in terms of the occupation 
of L fermionic modes, 



that it can be obtain from the 2L x 2L correlation matrix of some major ana 
operators. So that, the numerical cost only grows polynomially with the 
number of sites in the reduced system and the entropy is just the sum of the 
entropy for every mode, 




(2.42) 




(2.43) 
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2.3. Analysis of entanglement behavior. 



The results were obtained using numerical techniques to diagonalize the 
correlation matrix eq.(2.40). Due to the points just mentioned they 

coincide with the analytical treatment employed by different authors after 
our work. Jin and Korepin[JK04, IJK05] employed Toeplitz matrix proper- 
ties to diagonalize the correlation matrix and to get the thermodynamical 
behavior of the entropy. Peschel[Pes04a] got the same results using cor- 
ner transfer techniques to solve the problem. Keating and Mezzadri[KM04] 
applying random matrix theory also obtained the same conclusions. 

One of the main results in this work is the complete description of the 
behavior of a theory by the scaling properties of the entropy of entanglement. 
The fig. 2.9 shows this fact for the parameter space of the XY model. Points 
where the entropy diverges correspond to critical regions while the normal 
phase is characterized by a finite entropy of entanglement. Note and compare 
the general structure of this plot with the fig. 2.3. In principle, this figure 
appears quite qualitative but scanning through the different regions we will 
realize about the amount of information in this plot. 

Figure 2.9: Entropy of entanglement in a XY spin chain for L = 100 
contiguous sites as a function of the anisotropy 7 and magnetic field A. The 
regions where the entropy diverges correspond to the critical point of the 
theory. 



The line that corresponds to the XX model, i.e. zero anisotropy (7 = 0), 
is described by a gapless theory of two noninteracting free fermions. In the 
fig. 2. 10, we see that the entropy as a function of the number of sites L of the 
subsystem grows unbounded like a logarithm. The leading behavior does 
not change in the interval of the magnetic field |A| < 1 but some logarithmic 
correction are introduced when the field is turned on. In the limit case 
when A — ► 1 the entropy vanishes. We have seen in the diagonalization of 
this model that the magnetic field is a marginal quantity that just redefined 
the Fermi points of the theory, that is the reason why the entropy is only 
sensitive to logarithmic corrections in the value of the magnetic field. From 
the plot 2.10, the entropy is described by the relation, 

S (L, 7 = 0, A) = i log 2 L + i log 2 (1 - A 2 ) . (2.44) 

The error in the coefficients of the fit can be got with usual technique giving 
a value of the order 10 -8 . 

There is another critical region characterized by a free majorana fermion, 
that corresponds to values of the magnetic field A = 1. This region is de- 
scribed by the upper curve in the fig. 2.11 and the fig. 2.12 represents the 
correction with the anisotropy 7. So, in this region the entropy grows un- 
bounded. Its increasing and the leading order follows a logarithmic function. 
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Figure 2.10: Entropy for a reduced system of L spin in the XX model. The 
upper curve corresponds to the system without magnetic field. In the other 
curves, it can be seen that while the magnetic field increases, the entropy 
decreases. At the point of the field A = 1 the ground state is a ferromagnet 
and the entropy goes to zero. 



When we solve the problem, we saw that the anisotropy redefined the veloc- 
ity in the theory, so is a marginal correction, and therefore the anisotropy 
just induced logarithmic corrections. From these plots, the entropy is de- 
scribed by, 

5(L, 7 ,A = 1) = ilog 2 L + ilog 27 , (2.45) 

D O 

where the errors in the fit are of the order 1CP 8 . 

Away from the critical points, the entropy grows like a logarithm un- 
til the correlation length defined by the inverse of the gap in the energy 
spectrum, from that point the entropy is kept constant with the value, 

S(L,X^1) = ilog 2 -, VL>l/?n. (2.46) 



Figure 2.11: Entropy for the reduced density matrix of L spins in the Ising 
line, i.e. anisotropy 7 = 1. In this region, there is one fixed point at A = 1 
where the entropy grows unbounded with length of the subsystem. For mod- 
els with magnetic field A ^ 1, the entropy grows with a logarithm behavior 
till the correlation length, defined by the gap in the energy spectrum, at that 
point the curve is saturated by the gap. In the limit case A = the ground 
state of the system, if there is no symmetry breaking field, is described by 
a cat state, that it is the reason for the value one in the entropy. 



Figure 2.12: Logarithmic correction to the entropy in the XY critical line 
due to the anisotropy. The points fit the function | log 2 7 with errors of the 
order 10~ 8 . 



The plateaus in the entropy away from the critical point are studied in 
detail in the last two figures in the particular case of the Ising line, when 
7 = 1 with the variation of the magnetic field. In fig. 2. 13, a divergent point 
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at A = 1 shows the critical Ising point. From this point, there are two 
different behaviors in the plot. If the magnetic field increases the ground 
state of the system goes to a ferromagnetic state where every spin is aligned 
in the field direction and the entropy vanishes. If the magnetic field decreases 
to zero the theory is two-fold degenerate. If there is no symmetry breaking 
field, the ground state is the symmetric superposition of both vacuums, a 
cat state, and the entropy goes to one. But this solution is highly unstable, 
any small perturbation in the hamiltonian breaks this symmetry and chooses 
one of the vacuum. This scenario is described by the points in the lower part 
of fig. 2.13 which were obtained by DMRG (Density matrix renormalization 
group) techniques [Whi92] because introducing any symmetry breaking field 
make the system non integrable. 

Figure 2.13: Entropy of a reduced density matrix for L = 100 contiguous 
sites in an Ising chain with magnetic field A. In the left part of the figure, 
the upper points corresponds to the model without any symmetry breaking 
field. The lower points shows the entropy for a model that includes an small 
field that forces the symmetry breaking. 



In fig. 2.14, it is shown how the entropy scales in the Ising model as the 
magnetic field approaches its critical value. The plot shows that the entropy 
scales as, 

AS /sm «~-log,(l-A 2 ) (2.47) 
6 



Figure 2.14: Scaling of the entropy in the Ising model as the magnetic 
field approaches its critical values A = 1. The asymptotes correspond to 
AS Isin 9 ~ ±log 2 (1 - A 2 ). 



These figures, fig. 2.13 and fig. 2.14, will become an important part in 
the description of entanglement loss and renormalization group flows because 
as we will see in the next section, the flow followed by the Ising model while 
the magnetic field increases corresponds to a renormalization group flow. 

2.4 Connection with Conformal Field Theory 

Up to now, we have seen that the leading behavior of the entropy in a critical 
spin model is given by: Sl ~ k log 2 L, where k is a constant that does not 
depend on every detail of the hamiltonian but it seems to remain fixed for 
a given critical region. An answer of why the entropy has this behavior is 
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given by the structure of any one dimensional local theory, gapless and with 
a linear dispersion relation. Conformal field theory, (see [BZJ89, DFMS97] 
and references therein) is the area of research that studies the structure 
of these kind of theories and its physical consequences. Based in general 
arguments from conformal theory Holzhey, Larsen and Wilczek[HLW94] in 
1994 showed that the leading order behavior for the entropy in a bipartite 
scenario is given by: Sl ~ ^^log 2 L, where ^£ is the universal coefficient 
(called central charge) that we found for the spin chain, | in the XX and 
Heisenberg model and g in the XY model. Similar results were obtained 
by Korepin in 2004[Kor04]. Calabrese and Cardy[CC04] in 2004 extended 
the analysis to massive theories and different partition of the system. Also, 
scaling entropy has been studied in random quantum critical points by Refael 
and Moore [RM04] getting non-rational effective central charges. 

Conformal field theory has classified the properties of any critical field 
theory with linear energy dispersion in one dimension just from general 
arguments like the symmetry of the order parameter at the phase transition, 
leading to the concept of universality of the model, i.e. different physical 
realizations of a system can be described by the same theory. The coefficient 
c that appears in the constant term of the entropy is an omnipresent quantity 
in the classification of conformal field theories, it is the central charge or 
conformal anomaly. The central charge is one of the quantities that defines 
the universality and properties of a given critical system[BCN86, Aff86]. 
c = 1/2 corresponds to a free fermion theory and labels the Ising universality 
class which agrees with our calculations, c = 1 corresponds to a free boson 
theory. In fact, the continuum limit of the XX and Heisenberg model can 
be seen as a Sine-Gordon model, that is the theory of a bosonic field. 

Another remarkable meaning of the central charge is that can be seen as 
measure of the number of species of a given field in the critical region, i.e 
it is additive with the number of fields in a theory. So, a theory with two 
free fermionic fields has total central charge equal to one, this is the reason 
why c in the XX model is twice the one of the Ising case. This meaning 
links with the fact that the central charge appears explicitly in the equation 
for the entropy that is a measure of the number of degrees of freedom that 
describe the system. 

2.5 Scaling of entanglement in the Lipkin-Meshkov-| 
Glick model 

In this section, we will describe[LORV04] a model which is placed in a 
completed connected graph, that is, a model in which every site interacts 
with each other. This system allows us to study the role of the connectivity 
in the entanglement properties. For this purpose, we will study a spin system 
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2.5. Scaling of entanglement in the Lipkin-Meshkov-Glick model 



characterized by the Lipkin-Meshkov-Glick model[LMG65] (LMG) 

H = ~ Jj E {fi<fi + 7^) -h"£of (2.48) 

i<j i 

where cr" is the usual Pauli matrix at the j position in the a direction, 
N is the total number of spins in the systems. This hamiltonian describes 
a ferromagnetic model with an anisotropy in the xy plane given by the 
parameter 7 and a magnetic field h in the z direction. The LMG hamiltonian 
can be rewritten in terms of the total spin operators S a = \ Ylj a f > 

H = Hf L { {SZ)2 + y " ^ 2 ) " 2hSZ + ^ ((5+)2 + (s " )2) ■ (2-49) 

From this expression, it is straightforward to realize that the ground state of 
the system lives in the maximum spin sector S = y. Also, due to the high 
degree of symmetry of this model, a convenient basis to characterize the 
system is given by the Dicke states \^,m) that are fully symmetric under 
any permutation, they are eigenstates of S 2 with eigenvalue f (f + l) and 
S z with eigenvalue m. 

As we have just mentioned, the ground state of the system for the dif- 
ferent values 7 and h belongs to the permutational symmetric subspace. 
Therefore, the size of the Hilbert space grows in a polynomial way with the 
number of sites of the system. In fact, the reduced density matrix of L sites, 
from the ground state of ./V spins, is spanned by the set of L + l Dicke states, 
then the entropy of entanglement is always saturated by Sl < log 2 L + 1. 
This properties make that this model can be treated efficiently with numer- 
ical methods. 

In the fig. 2.15, we show the general structure of the ground state char- 
acterized by the von Neumann entropy of a bipartition of the system. As 
in the case of the XY model, we will scan and analyze every region of the 
parameter space in the following. 

Figure 2.15: Entanglement entropy for N=500 and L=125 as a function of 
h and 7. 



In the isotropic limit 7 — ► 1, the LMG model is diagonal in the Dicke 
basis. The ground state for a magnetic field h is given by the Dicke state 
|f ,M), its energy is given by £? (^7 = 1) = -f + M 2 ^ - 2hM and 
magnetization in the z direction, 
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where I[x] is the round value of x. 

If instead of the magnetization M, we use the number n = M + ^ of 
"up" spins to describe the ground state, i.e. |y,M) = \N,n), we can see 
that any bipartition of the system is quite simple to realize, 

N 

\N,n)=^2^pi\L,l)®\N-L,n-l), (2.51) 
;=o 

where the subsystems are blocks of size L and N — L sites and pi is the 
hypergeometric distribution, 

'L\ (N -V s 

= 7 ]v^ ' ^ ' 

n 

So, the entropy of any subsystem of L sites is given by 

L 

SL,N(h,7) = - ^2pilog 2 Pi- (2.53) 

1=0 

In the limit iV > L > 1, the distribution p; approaches a gaussian one 
;h mean va 
is recast into, 



with mean value I = nj^ and variance a 2 = n ^ N N n ' N N L j^, and the entropy 




W, = .)H r "- " / + (I"*") (2 .54) 

rt > 1 

where in the interval h G [0, 1) we keep the sub leading term N ^ L to make 
explicit the symmetry Sl,n = Sn~l,n and in the region h > 1 the ground 
state is fully polarized in the direction of the magnetic field, so it is separable 
and its entropy is zero. 

In the anisotropic case, i.e. 7 7^ 1, the LMG hamiltonian is not diagonal 
but still can be expressed in terms of the Dicke basis, and as we have seen, 
the cost to get the coefficient of the ground state is just polynomial with the 
number of sites. 

Figure 2.16: Entanglement entropy at 7 = as a function of h for different 
values of N and L. Outside of the critical region, the entropy only depends 
on the ratio jj. 



In the fig. 2.16 and 2.17, we plot the dependence of the entropy with 
magnetic field. Several results appear in these plots that are worth to point 
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out. From the first figure, it can be seen that the value of the entropy 
for a given value of h and 7 only depends on the ratio -fe. In the limit 
h S> 1 the ground state of the system is totally polarized in the direction of 
the magnetic field, and then, the entropy goes to zero. When h — ► 0, the 
system is in the symmetric superposition of the polarized vacuums in the 
±x direction, i.e., the ground state is described by a "cat" state or GHZ- 
like[GHZ89] state, so the entropy goes to the value one. As the magnetic 
field approaches the critical value, h — > 1, the entropy displays a logarithm 
divergence 

S Lj jV~-ilog 2 |l-fc| (2.55) 



Figure 2.17: Entanglement entropy as a function of h near the critical 
point for 7 = 0. The full line corresponds to the fitting law S^^/i, 7) ~ 
-|log 2 \l-h\. 



In the plot 2.18, we display the scaling of the entropy with the number 
of sites L that are considered in the reduced density matrix. The law that 
fit the curves is given by 

S L>N (h = 1, 7 + 1) ~ i log 2 U^j^) ( 2 - 56 ) 

with subleading corrections given by the anisotropy parameter 7. These 
corrections are studied in the next figure 2.19, from which it can be seen 
that they are logarithm correction given by 

S L A h = 1,7)- S L>N (h = 1, 7 = 0) ~ \ log 2 (1 - 7) (2.57) 



Figure 2.18: Entanglement entropy as a function of L at the critical point 
for different 7 and N = 2000. The full line corresponds to the fitting law 

52 AT 



Figure 2.19: Entanglement entropy at the critical point h = 1 as a function 
of 7. The full line corresponds to the fitting law SL,N{h = 1,7) — SL,N(h = 
l,7 = 0)~±log 2 (l- 7 ). 
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This detail analysis of the entanglement entropy of the ground state in 
the LMG model shows an incredible resemblance with the XY model in one 
dimensions, although the former is placed in a completely connected graph 
and the latter in a one dimensional lattice. Due to the symmetry properties 
of the LMG, we have argued that the entropy can only scale with a logarithm 
behavior, but it seems that a deeper understanding of scaling in the LMG 
model is still needed. 

2.6 Conclusions and related works. 

2.6.1 Critical and non critical spin chains. 

The core of this work analyses scaling properties of the entropy of entangle- 
ment of the ground state in relevant spin systems. With this detailed study, 
we have been able to describe the phase diagram of these systems. This 
analysis and the characterization of the different phases can be linked with 
concepts like symmetry and universality that appear at the critical points. 
So, the main results of this sections are summarized as follows: 

1. There exists universal scaling law of entanglement in 1+1 dimensional 
systems in the critical regions. 

2. Entanglement behavior is controlled by conformal symmetry. 

3. Away from the critical points entanglement is saturated by mass scales. 

As a consequence of these three points and understanding how the size of 
the Hilbert space of a subsystem scales, it is straightforward to understand 
the success of DMRG (Density Matrix Renormalization Group) or similar 
variational methods to describe and implement the ground state of one di- 
mensional systems. In this line, a main question is to get an efficient way to 
simulate quantum systems in higher dimensions [VC 04]. 

2.6.2 Entanglement in higher dimensions. 

Table 2.1: Scaling of entanglement entropy of a subsystem with a typical 
length R. 



Entropy Sr 


Critical region 


Non-critical region 


d = 1 dimension 


log 2 R 


log 2 m 


d > 1 dimensions 


R d -' 



One of the first motivations to study entanglement properties in quantum 
systems was to understand better the origin of black hole entropy [Bek94] 
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and the area law, i.e. in a three dimensional system, the entropy of a 
region, with a typical length R, is expected to scale as S oc R 2 , in a d- 
dimensional system, S oc R d ~ l . In this context, many works have try to link 
this entropy with the entanglement entropy between the inside and outside 
mode of the event horizon (see for example [BKLS86, CW94, KS94, FPST94, 
Kab95]). Srednicki in 1993[Sre93] using a Klein-Gordon model, explained 
the area law behavior of the entropy of entanglement as a consequence of 
the local properties of the interactions and the Schimdt decomposition of 
the ground state of the system between inside and ouside modes of some 
region, i.e. <p = Yl I 1 4>in (pout- Because the Schimdt parameters are the 
same for both subsystems, the entropy of entanglement is the same seeing 
from inside or outside. So, if there is any dependence in some geometric 
property of the subsystem, it should be the boundary that they share. Prom 
a quantum information perspective, better bounds and analytic studies have 
been developed for bosonic[PEDC05] and fermionic fields [W6105], arriving 
at similar conclusions. 

2.6.3 Non-local interactions and dynamical evolution of en- 
tanglement. 

The study of entanglement in non-local systems has been used to describe 
how entanglement appears in quantum gases or spin glasses [CHDB05]. This 
kind of systems gives a model to characterize decoherence. Even, pos- 
sible consequences to adiabatic quantum computation success arise from 
these analysis [LO04, OL04]. There, it has been shown that adiabatic quan- 
tum computation uses exponentially many entanglement when solving NP- 
complete problems, which has direct implications in classical simulation pro- 
tocols. Also, it has been applied in the context of dynamical evolution of 
quantum computers with imperfection [MBF03]. Finally, it was noted that 
the dynamical evolution of entanglement is constrained by properties and 
symmetries like causality [CC05] 



Chapter 3 

Entanglement and RG-flows 



The renormalization group (RG) is a method designed to describe how 
the dynamic of some system changes when we change the scale (distance, 
energy,...) at which we probe it. A theory at a larger scale employs only a 
finite part of degrees of freedom from theories at smaller scales, as the irrel- 
evant details are integrated. The decoupling of physics at larger scales from 
the one a smaller scales is the reason why it is possible to do atomic physics 
without knowing the final structure of the constituent of the nuclei. So, 
the description of low-energy physics should not need a detailed knowledge 
of the laws of Nature at very high momenta. The transition of informa- 
tion from the smaller scale theories to larger scale theories is irreversible, 
since in such transition we integrate out many irrelevant degrees of freedom. 
The general aim of renormalization group method is to explain how the 
decoupling between different scales takes place and how the information is 
transmitted from scale to scale. 

A concrete scheme for renormalization was shown by K. Wilson[WK74, 
Wil75]. One of the main objects in his picture was the space of theories S, 
as the space of all possible hamiltonians for a given field. So, let /xi, /i2 be 
some scales of momenta or energy, with ^2 < A*i- For each theory Ti £ S 
there is another theory, TZ^ lfJj2 7i € S which is the large scale theory (or 
effective theory), at the scale fj,%, for the original one "H at n\. Thus, there 
is a map 72. Ml/i2 : S — ► 5, with ^ 1 ^2^M2M3 = ^1/^3 that defined the action 
of a semigroup usually called renormalization group. 

The aim of this second work is to analyze how correlations between 
short and long distance degrees of freedom are lost in RG transformations. 
To get this goal, we employ tools from Quantum Information to measure the 
degree of entanglement in the vacuum of a theory. Mainly, we will perform 
a quantitative analysis in quantum Ising model with a transverse magnetic 
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3.1. RG transformation of free theories. 



field as the paradigm of a gaussian fermionic system. In the last part, we 
will describe the flow of the theory in the space of hamiltonians just with 
quantities like entropy or majorization. 

3.1 RG transformation of free theories. 

Any renormalization group transformation [Sha94, Mor94, BTW02, BHLM95]| 
has two steps: 

1. Integration of high energy modes that let a theory that describes the 
same long-distance properties. 

2. Rescaling of parameters that allows to compare the initial with the 
renormalized theory because at this point they share the same set of 
degrees of freedom. 

The action of the renormalization group is particularly simple in free theories 
like the XY model. We have seen that the XY model can be seen as a free 
theory of fermions, 

H = # Vfij, ~ [ 1/a d<f> Jm 2 + ^ 2 +aV b\b^ (3.1) 

J—l/a J — 1/a 

where we have explicitly regularized the theory with a lattice spacing a 
or ultraviolet cutoff 1/a, and in the second equality, we have developed 
the dispersion relation to order o(t/> 4 ). Due to the fact that this theory is 
noninteractive, the high energy modes are decoupled from the low energy 
ones, i.e. the partition function can be split in two independent factors: 
Z = lim0_ +oo ire~~P H = Z\ k \ <e -i i a Zu.\ >e -i / a , where I is the scaling parameter 
and the hamiltonian for the low energy modes is written as, 

%|< e -7a = f ^ # Jm* + v 2 F cp + a*cl>* b\b^. (3.2) 

At this stage, it can be seen that in these kind of theories all the non- 
trivial properties in the RG transformation come from the rescaling of the 
parameters that is implemented as follows, 

4> = e- l <j)\ H = e~ l KH, b 4> = e xl b 4>l . (3.3) 

The first two equations take into account the rescaling of the momenta and 
energy as we change the scale we are using, while the second takes into 
account the rescaling of the field, where the scaling dimensions x are fixed 
by the invariance of the kinetic term vp(j) 2 , in this case x = ^. Finally, the 
renormalized hamiltonian is recast into, 
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So, the theory can be seen as the initial one where the parameters of the 
hamiltonian have been renormalized, i.e. m — > e l m and a — > e~ l a. Then, 
an infinitesimal transformation in the scale, implies 

8vf dm da 

These three equations describe how the parameters of the theory change 
under the renormalization group, they are usually called beta functions. 
Successive RG transformations define a RG flow in coupling space. In this 
example, we see the three kinds of behavior in the parameter space of any 
theory: 

1. Relevant parameter.- An infinitesimal perturbation of the parameter 
grows under successive RG transformations. So that, the perturbation 
is important or relevant in the long distance behavior. In our example, 
it will correspond to the mass term which drives the system from the 
Ising fixed point towards an infinitely massive theory. 

2. Marginal parameter.- A perturbation in the parameter does not change 
under RG transformations. In our example the kinetic term. 

3. Irrelevant parameter.- A perturbation in the parameter will decrease 
under RG. So, it is irrelevant for the physics at long distance and if 
we are looking for an effective theory, we can neglect its effect. In our 
model, this is the case for the parameter a. 

A fixed point in an RG transformation is defined where the parameters 
of the theory remain constant. In free theories, the fixed point appears when 
a and the mass goes to zero, i.e. at the critical points. 



3.2 Entanglement loss and RG-flows. 

In one dimensional systems, entanglement loss in RG-flows [LLRV05] can be 
analyzed in three different steps that we described in what follows. The 
first step is just a direct consequence of the structure of critical theories 
of one dimensional systems and so, it is quite general. For the second and 
third steps, we will fix our attention in massive fermionic theories to develop 
all these ideas, taking the Ising model as a representative of these kind of 
models. 

Global loss of entanglement.- As we mentioned, conformal field theory 
describes the properties of one dimensional critical systems and looks for its 
consequences. One of these consequence was shown by Zamolodchikov[Zam86]| 
in 1984, when he stated the C-theorem (see also [CFL91, FL98a, FL98b, 
AnsOO]). There, he proves that unitary and Poincare invariant theories that 
are perturbed with a relevant deformation from a critical point, under RG 
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transformations, they flow to another fixed point. This flow is characterized 
by a function C(gi,/j,) called c- function, that depends on the parameters of 
the theory g^ and the scale at which we test it, fi. This function is monoton- 
ically non-increasing, i.e. — (u^J ~5g~C — an d remains constant at the 
fixed points of the renormalization group transformations. The fixed value 
of the function built by Zamolodchikov coincides with the central charge 
of model at the critical point, i.e. C(g*,fi) = c. Therefore, the central 
charge in any ultraviolet fixed point, i.e. when all the details, relevant or 
irrelevant, are still present, is greater or equal than the central charge of 
the infrared theory or effective theory, i.e. c uv > c IR . In the last sec- 
tion, we showed that the central charge of a critical point is directly related 
to the behavior of the entropy of entanglement in the ground state of the 
theory, S(L,g*) = log 2 L. Then, the entanglement entropy in the UV 
theory is never smaller than the one in the IR theory, i.e. S uv > S IR . 
So, the first link between RG-flows and entanglement appears as a conse- 
quence of the structure of conformal theories and relates the entropy between 
two fixed points. This fact tries to explain the C-theorem by an entropic 
reasoning[FPST94, CH04] and points at the loss of entanglement as the 
cause for the irreversibility in the RG flow. 

Monotonic loss of entanglement. - Looking for a deeper description of the 
loss of entanglement under RG transformation, at this second step we will 
follow the RG trajectories in the Ising model and check that this loss appears 
step by step in the flow. Our analysis comes from the observation (fig. 2.13) 
that entanglement entropy itself behaves as a c-function in free fermionic the- 
ories, so it is a monotonically non- increasing function in the flow of the mass. 
In fact, we saw that close to the critical point, AS(L, A) ~ — g log 2 (l — A 2 ). 
From the plot we observe that there are to limit points, if the magnetic field 
A >> 1 then the entropy goes to zero but in the phase where A — ► the en- 
tropy approaches the unit. This value is characteristic of cat states of GHZ 
states which are a symmetric superposition of two orthogonal and polarized 
ground states, i.e. = 4j (| |T • • • T) + I II • • • [))■ Nevertheless this state 
is highly unstable and any small perturbation of the system, makes it fall in 
one of the orthogonal vacuums. Then, if an small symmetry breaking field 
is applied in the hamiltonian, the entropy will approach to zero, where the 
ground state is totally polarized and separable. 

Fine-grained loss of entanglement.- (see also [Oru05]) This third step 
shows that, at least, for fermionic gaussian theories, the loss of entanglement 
entropy under RG transformations is just the consequence of the fact that 
initial and final states of the system are linked by majorization relations (see 
appendix A). As we will sketch, these transformations are deeply connected 
with unitary concepts which links with the fact that in the C-theorem one 
of the main assumptions is the unitary of a theory. 

In the second chapter, we have seen that the reduced density matrix of 
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the Ising model can be written in terms of normal fermionic modes, i.e. 

p= n p*- ( 3 - 6 ) 

is the reduce density matrix for the <j> mode, that written in terms of the 
occupation numbers, is recast into 

^ = ^(|0}(0|+e-^W|l)(l|), (3.7) 

where Z^ is a normalization factor, uj^ is the dispersion relation for the 
4> mode (their explicit expression for the XY and XXZ model was given 
by Peschel and collaborators [PKL99, Pes04a]), A is the magnetic field that 
characterizes the behavior in the Ising model. 

In what follows, we will see which are the conditions in the dispersion 
relation in order to fulfill majorization relations (see appendix A) mode by 
mode and at every infinitesimal step of the RG flow. 

The density matrix majorizes p^, written p^ ~< p^, if and only if they 
are related by a positive, unital and trace-preserving map, i.e. a double 
stochastic transformation. In terms of the vector of eigenvalues of each 
density matrix implies that the vector A(p^) is given by a probabilistic 
combination of permutations of the vector of eigenvalues A(p^), i.e. 

A(p*) = Y,pf P ^Pt>) = D<A(Pt), (3-8) 



where is a double stochastic matrix, the real numbers p^ are such that 

< p^ < 1, ^2jP& = 1 and {Pj} is a set of permutation matrices. In the 
case of fermionic modes, the vector of eigenvalues has just two entries and 
there are two permutation matrices, 



o 



(3.10) 



then, any two fermionic density matrix are linked by majorization relations 
iff 

J_ / 1 \ _L_ / 1 e- Q 

{?-"*) ~ V e "" 1 

these equations are equivalent to 

(o) _1 + e"^ 1 - e -^e"^ 




<t> i _|_ e _a ^ 1 — e -2 ^ 
(i) _ 1 + e"^ e"^ - e" <i 
<t> ~ i -\- e _tJ 1 — e" 2 ^ 



(3.H) 
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with the constraints py +p^ = 1, < p? < 1 and < p^ < 1. 

Parametrizing any flow in the space of the density matrix by the variable 
r, the dispersion relation will be a function of this parameter uu = ujAt). 
An infinitesimal transformation that brings r — > f = r + Ar, will change 
the dispersion relation to = Co^ — ^^-At + O ((At) 2 ) , where u)^ = u^(f ). 

Then, a simple analysis shows that majorization holds mode by mode at 
any infinitesimal step in the flow of r iff 

dCj<h At 

0< — ^<1 3.12 

- Qt - e'^ ~ 

from this equation, it can be seen the three aspect in the flow and the 
dispersion relation that have to be checked: (?) the increasing behavior of the 
dispersion relation from ^§jr', (ii) the monotonic behavior of the parameter 
t; {Hi) the sign of the dispersion relation. 

Usually, the dispersion relation is positive defined, so that, the main 
contain, for ^ lm plies that, sgn(^) = s g „(A T ) 

Finally, when majorization holds mode by mode, the whole density ma- 
trix of the system, also fulfills these relations. This statement is due to 
the fact that the tensor product of double stochastic matrices defines a new 
double stochastic matrix V, 



A(p) = JjA(p^) = n^ A (^)= \ UdAa(p)=VA(p). (3.13) 




As an explicit example, we can see that the Ising model fulfills every 
conditions in the RG flow in order to relate the eigenvalues of the density 
matrix with majorization relations and so the loss of entanglement is the 
consequence of this fact. 

We have seen that the RG flow in this model is characterized by the 
increasing of the mass term, Am > 0, so this variable parametrizes the flow 
of the model, i.e. r = m. We also know that in the Ising case, m = |1 — A|, 
where A is the magnetic field. Then, the departure from the critical point can 
be seen as the change of the magnetic field around A = 1. Peschel, Kaulke 
and Legeza[PKL99] showed that the dispersion relation in this model is: 



j (yi-(i/A)2) 

(3.14) 



(2^ + 1)tt vv /(1/a) ; A>1 



2<jm VV J(A) ' A < 1 



where I(x) denotes the complete elliptic integral of the first kind[GR94]. 

From the explicit expressions of the dispersion relation for the Ising 
model, it is quite simple to see that —g^- > and so majorization hold 
for every mode eft. There is just one last comment about the zero mode 
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Figure 3.1: Dispersion relation uj(\) for the Ising model. The first plot shows 
the function when A < 1 while the second plot is for A > 1. 



in the unbroken phase, when A < 1. The eigenvalues of this mode are 
A(po) = \ (1> 1) independently of the value of the magnetic field. This 
special behavior produces a cat state in the limit A — > 0, i.e. in a local 
spin basis = A= (| || • • • |) + | J.J, ■ ■ ■ j)). This state is quite unstable to 
any perturbation in the hamiltonian and in fact it violates the clustering 
principle. If a symmetry breaking field is applied to the hamiltonian, the 
state will flow to one of the polarized vacuum. 

3.3 Conclusions and outlook. 

We have seen that a deeper analysis using tools from quantum information 
to describe and characterize pure quantum correlations has revealed a highly 
ordered structure in the vacuum of the studied models. This structure seems 
to imply an order relation in the space of states under scale transformations 
which we use to point at the entanglement loss as the root of the irreversibil- 
ity in the renormalization group flows. So, two main ideas appears from this 
study: 

1. Entanglement is non increasing under RG flows. 

2. There exists a close relation between majorization and C-functions. 

At least two obvious points remain as open questions, the generalization of 
these links to interactive theories and if majorization can give new insight 
to higher dimensional systems. 



Chapter 4 

Matrix Product States. 



In the last decades, one dimensional models have attracted much at- 
tention and effort to understand their behavior. Several reasons motivate 
this effort: many physical systems can be mapped to these models and new 
collective phenomena appear in low dimensional systems. In this context, 
a new sets of methods have been developed to simulate in an efficient way 
these models and characterize their ground state. The methods that we 
describe in what follow, share as a building block a set of local finite ma- 
trices that implement the information of the system. We will show how to 
build these matrices, their properties and how can be used to compute, in 
an efficient way, any correlator. At the end, we will see how to perform 
scale transformations within these formalism [VCL + 05], how to realize an 
exact real space renormalization group at the level of the ground state of 
the system and how it can be found the fixed points of the transformation. 



4.1 Definition and properties. 

To define the matrix product state, we will show two possible ways to built 
them. The first one is in the spirit of the Density Matrix Renormalization 
Group (DMRG) method of S. White[Whi92] (see also [Sch05]), and it is 
based in the Schmidt decomposition of a quantum state (see appendix A). 
The second one, the older point of view, comes from the fact that the ma- 
trices can be seen as a bond between the neighbor sites in a chain. This one 
was the idea of the valence bond model of I. Affleck, T. Kennedy, E.H. Lieb 
and H. Tasaki[AKLT87, AKLT88], generalized by M. Fannes, B. Nachter- 
gaele and R. F. Werner [FNW92]. The link between both point of view was 
initially given by S. Ostlund and S. Rommer[OR95]. The results from scal- 
ing of entanglement allow us to explain the success of these formalism in 
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one dimensional system. Finally G. Vidal[Vid03, Vid04], F. Verstraete, J.I. 
Cirac, J.J. Garcia-Ripoll and D. Porras[VGRC04, VPC04b] have shown, in 
different works, how to generalize these methods to nonlocal interactions, 
different boundary conditions and how to perform dynamical simulations of 
quantum systems with classical resources. 

Following the fig. 4.1, the Schmidt decomposition of a quantum state 
between the first site and rest of the chain reads, 

l^)=E^|aS 1) )|aS 2 "- JV) ) (4.1) 

ai=l 

where is an orthonormal basis of the first subsystem and {|a^ 2 , " Ar - 1 )} 

is an orthonormal set for the rest of the chain, the real numbers | < 
Men < 1) Sai^ai) 2 = 1} are the Schmidt coefficients. Finally, the integer 
D is the Schmidt rank in the decomposition that, for easy writing, we fix to 
the maximum of any bipartition. If we express the state of the first site in 
some initial local basis, then the state is written, 

l*)=EE IsiK^rf) \«?- N) ) = E E \<$~*) (4.2) 

ai=lsi=l cki=1si=1 

where in this case, d is the dimension of the local space and = /j, ai (si\a^)\ 
is the tensor for the first site. 

Figure 4.1: Schematic first step in the construction of the matrix product 
state using Schimdt decomposition in a chain with iV sites. 



Carrying the process to the second site, 

D d 

l*> = E E \ s 2)A s a \ a2 \^- N) ), (4-3) 

{a}=l{s}=l 

where the sums {a} and {s} are made over the different configurations of 

(2) 

a>i, ct2, s\ and S2 and the matrix A s * ia2 = fi a2 {s2\aiQ. 2 )■ Then for the 
whole chain, the state is written as follows, 

D d 

l*>=E E k^'H^CVJ^- (4-4) 

W=l{s}=l 

Looking at the bulk of an infinite chain, the basis in the bipartition for 
the right subsystem will be written as, 

D d 

I«l-i)= E E ^A-vl-^W), (4.5) 

<*£, = ! s L = l 
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where \a L ) = Eq l+1 =i Es l+1 =i \sL+i)A s a L L + a 1 L+1 \a L+1 ) and if {|a L _i)} is an 
ortho-normal basis then 

d 

(A 3 ^ A SI < =1, (4.6) 

SL = l 

which defines the condition in the A s matrices to be a positive, trace pre- 
serving map. Within this approach, the matrices A SL represent the change 
in the description of a system when it is considered one site more, they act 
as a transfer matrix. Every time we insert one site more in the subsystem, 
a new matrix A SL link it with the rest of the chain and, in that way, a suc- 
cession of matrices {^4 Si } is produced. If this succession has a limit, then 
we could represent the final state of the chain by this set of matrices. 

So, this method is computationally efficient if the Schmidt rank D is 
bounded. But this situation appears usually in one dimensional system 
with local interactions. When we perform the bipartition of the system 
between L sites and the rest of the chain, the Schmidt rank measures the 
number of degree of freedom needed to describe the subsystem and is always 
bounded by the entropy of entanglement, i.e. D > 2 Sl . Using the analysis 
of the scaling of entanglement in one dimensional models, we know that for 
finite systems of length N the maximum of the entropy goes as Sl ~ log 2 N 
and for infinite gapped systems, it goes like Sl ~ log 2 — . So, in any of 
these cases the Schmidt rank is bounded by a constant. Even at the critical 
points, we know that the entropy grows with the size of the subsystem L as 
Sl ~ log 2 L, then, although the Schimdt rank is not saturated in this case, 
its growth is polynomial with the size of the subsystem, D ~ L. 

In a translationally invariant state, the whole set of matrices can be 
represented by just one matrix independently of its position, i.e. {A SL } ^ 
A s . Then, a translationally invariant state can be described by the product 
of this matrix as follows, 

d 

|¥) = ^2 tr (5 ■ A Sl ■ ■ ■ A SN ) \s u ■■■ ,s n ) (4.7) 

{s}=l 

where the matrix 5 implements the boundary conditions of the system. In 
particular, if the state has periodic boundary condition this matrix is the 
identity, 5 = 1. 

The second point of view, for the matrix product state, comes from 
the identification of the lower indexes of the matrix A s a p with two ancillary 
subsystems used to implement the state of the physical system s. So, this 
matrix appears as a projector from the Hilbert space TiS a ' of the ancillae to 
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Figure 4.2: Valence bond picture of a matrix product state. The first line 
with dark blue spots represents the physical state of the spin chain, while 
the second line represent the implementation of the state with two ancillae 
systems per site and a maximally entangled state between neighbor sites. 

the real one 7i iphys \ i.e. 

A:C D ®C D -^C d , 

* (4.8) 
\a,P) -+A\a,0) =2^A s ap \s). 

s=l 

In this way, the correlations in the chain (fig. 4. 2) are implemented by maxi- 
mally entangled states in the ancillae subsystem of neighbor site, i.e. Yla=i \ a ) I 
and the final description of the state is, 

D D D d 

|*)=Aj>)|a)£|/3>|/?>= I>Mitfl*>l0>- ( 49 ) 

a=l (3=1 o,/3=l s=l 

This one was the original method due to Affleck, Kennedy, Lieb and Tasaki 
(AKLT) to built the exact ground state of antiferromagnetic system of a 
spin one chain using the symmetric subspace of two spin 1/2 ancillae. 

Once the state of the system is built, it is important to know how to 
get any correlator from it. In the matrix product state formalism, this 
task is implemented defining a transfer matrix . For instance, to get the 
normalization in a translational invariant state, 

d d 

(*|*)= E (si-s N \tr(A Sl ..A SN )*tr(A il ..A SN ) \h~s N ), (4.10) 

M=l {8}=1 

defining the D 2 xD 2 transfer matrix E = Y7 s =i{ As )*® AS \ the 

normalization 

is just (^l^) = tr (E N ). If the state is normalized, then, the spectrum of 
the matrices E is such that {A^J |AJ < 1, 1 < H < D 2 }. If the state 
fulfills the clustering principle 1 , then, the maximum eigenvalue Ai = 1 is 
unique [FNW92]. 

In this representation the two point function of any two operators at the 
site (i) and (j) is written as 

(O«C>0')) = tr (e^+^O^E^-^O^A (4.11) 
lr rhis principle requires that lim^^^ (*|O(x)O(0)|1') = (*jC>(x)|*)(*|O(0)|*} 
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where O = X^{ s }=i S{s}=i(^ s )* ®A s (s\0\s). If we consider an infinite chain 

D 2 

and taking the diagonal form of E = |1)(1| + XL=2 •Vl/ x )(A t l> then, 



i) 



D 2 '— i— 1 



] (i) |l) 

(4.12) 



n=2 



where £^ = log ^ | defines the correlation length. From this expression, we 
see that the long range behavior of the system is implement in the eigenvalues 
of the transfer matrix E, that is invariant under the redefinition of the 
local basis. An important point that appears about the matrix product 
state formalism is that any correlation function is expressed as a sum of 
exponential. This implies that within this formalism to get a non exponential 
behavior in the correlation, i.e. in the scaling region or at distance between 
the lattice spacing a and the correlation length £, a < x < £, the number of 
eigenvalues of the matrix E should be as large as the number of coefficients 
in the Laplace transform of the polynomial decay in the correlation. 



4.2 Renormalization group on quantum states. 

In the last section, we have seen that matrix product states is a simple for- 
malism to implement translationally invariant state, that any expectation 
value can be calculated in an easy way and, also, they are optimal to min- 
imize the ground state energy. One of the reason for these properties is 
that the set of local A s matrices define a transfer matrix between neighbor 
sites. This fact allows to split the long and short range behavior, i.e. there 
is an explicit separation of scales. Then, the matrix product state are prop- 
erly suited to implement the classical Kadanoff[KGH + 67] block spin at the 
quantum level. 

4.2.1 Scale transformations 

Given a one dimensional quantum state |^), it can be decomposed, up to 
some local unitaries, by a basis {\s)} of the local Hilbert space and a set of 
local matrices {^4 S }. In this way, we define an equivalent relation ~ between 
two states ~ j^) iff they are related with some local unitary operation 
I^O = n®i=i Ui\i&). This equivalent relation is motivated by the fact that 
long range behavior does not depend on the definition of the local basis. 
As we have seen, the correlation length is defined by the eigenvalues of the 
matrix E = X^s=i(^ s )* <8) A s that is independent of the local basis. 
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A scale transformation in the matrix product state is denned by blocking 
neighbor sites and identifying the local Hilbert space as the space spanned 
by the states in the block, redefining the A s matrices and relabeling the 
position of the blocks: 

D min(d 2 ,D 2 ) 

^«.«+0 = ^A^A'ff- 1 = Yl (' ' : X Sj \ ',U (4.13) 

0=1 Sj=l 

where U and V are unitary matrices and {X s j } the coefficients of the singular 
value decomposition of the block spin A S2i A S2i+1 . The first equality in this 
equation is just the coarse- graining or block spin while the second part takes 
into account the rescaling of the sites in the chain and the identification of 
the relevant Hilbert space of the block. It is this second step the origin 
of the irreversibility in the transformation, because once it is identified the 
relevant space in the block, any information about the local sites is washed 
out. 

So, the renormalization transformation is implemented by: 

A s RGA S = X S V S (4.14) 

The dimension of the local Hilbert space in this transformation is always 
bounded by D 2 as it comes from the singular value decomposition of the 
product of two A matrices. 

As an example of this transformation, we will study how evolve the 
ground state of the AKLT model for a spin one chain. The AKLT model 
is defined by an antiferromagnet spin hamiltonian with Heisenberg-like in- 
teraction between neighbor sites. It is known that the ground state of this 
model is given by a valence bond state which is an example of the matrix 
product state. If we take as a local basis a combination of eigenvectors of 
the S z operator, such that, 

|*> = 4= (l+> + |-»; |y> = 4= (1+) - |-» ; \z) = |o> 

V2 ^V2 (4.15) 

S z \n) = n\n),n = {-1,0, +1}, 

then, the A s matrices are written in terms of the Pauli matrices, 

A x = ^-a x ; A y = -^a y ; A z = -^a z ; (4.16) 

from which the correlation length of this state is £ = In this case, 

the scale transformation is quite simple due to the properties of the Pauli 
matrices, o a a^ = 5 a pa° + iea^a 1 , where e a/ g 7 is the Levi-Civita tensor, a 
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is the identity matrix and {a, [3, 7} = {x, y, z}. So, 

A x A x \xx) + AyAy\yy) + A z A z \zz) = \a° (\xx) + \yy) + \zz)) 
A x Ay\xy) + AyA x \yx) = V (\xy) - \yx)) 
AyA z \yz) + A z Ay\zy) = V (\yz) - \zy)) 
A z A x \xy) + A x A z \yx) = % -a v (\zx) - \xz)) 



(4.17) 



3 



relabeling the Hilbert space in the block, 



(4.18) 



|6> = ^=(\xx) + \yy) + \zz)) 
\x) = -j= (\yz) - \zy)) 

\y) = -j= Qzx) - \xz)) 
\z) = (\xy) - \yx)) 
redefining the matrices A s , 

io = 1 ^ = n/2 lv = i ^l CT y A z = ^a z . (4.19) 
y/Z 3 3 3 K ' 

In this way, we have an exact representation of the state in a long scale, 
that, once we forget (project) the local spin state that define the Hilbert 
space of the block, the irreversibility in the transformation appears. In fact, 
the new local Hilbert space has physical meaning, they are real quantum 
state of the system: |0) is the singlet state and \x), \y), \z) form the triplet. 
At any point in the scale transformation the local state is written by, 



= y/i -3// 2 <t°|0> + ifi Yl aS l s ) ( 4 ' 2 °) 

s={x,y,z} 

with correlation length £ = ; og |i"-4^| • Then, a step in the RG flow is defined 
by 

1*04) - |*(£)>, (l-4/, 2 ) 2 = (l-4A 2 ) (4.21) 

As a consequence of this analysis, we can seen how the entropy evolves for 
any given number of sites L and correlation length £ = i og |i"-4^ 2 | ' 



S(L, n) =3 log 2 



2\2 L 



r (4-22) 
l + 3(l-4/, 2 ) 2 " i + 3 (l-4 M 2 ) 2i 
: !og 2 
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whenever L — > oo or £ — > 0, the entropy goes to S — > log 2 4, see also 
[FKR04],this fact is a consequence of the structure of the valence bond 
state produced by a topological hidden order [GA89] and that it is reflected 
in a non zero entropy, an infinite entanglement length[VPC04a, VMDC04, 
PVMDC05] or long range correlation in the string order parameter [dNR89]. 

4.2.2 Fixed points in the scale transformation 

Up to now, we have seen how to realize exact renormalization group trans- 
formation on quantum states. Usually, if limit cycles[GW02, MN03, LS04] 
are not considered, this kind of transformations end in a fixed point that 
characterizes the long range behavior of the state. So, determining the con- 
ditions for the fixed point in the transformation and analyzing the properties 
of the quantum state will allow to know better how the low energy properties 
emerges from a given state. 

In the last section, we based the transformation in the A s matrices but, 
as we saw, their explicit form depend on the choice of the local basis. Also, 
we found that the transfer matrix from these matrices E = y~] (A s )* <g> A s 
are basis independent. In fact, it is this transfer matrix what is used in the 
calculation of correlation functions and its eigenvalues determine the long 
range behavior of the state. The definition of the scale transformation with 
this operators is even simpler: E = E 2 , and so, a fixed point in this trans- 
formation implies studying the class of operators {-Eoo} = {h m »woo E n } or, 
in other words, analyzing the operators E such that the absolute value of its 
eigenvalues are zero or one. To carry this program, we will fixed our atten- 
tion in states that can be decomposed with matrices of dimension D < 2 and 
with the knowledge of completely positive trace-preserving maps[RSW02] we 
will give a complete classification of fixed points in these transformations. 

1. The first example of fixed point in these RG transformations is given 
by product states, such that A 1 = 1 and A s = Vs ^ 1. The trans- 
formation is given by A 1 = A 1 A 1 and the Hilbert space of the block 
is one dimensional |1) = |11). Obviously, E = E. 

2. GHZ states [GHZ89] or cat states is another example of a fixed point. 
The ground state of the Ising model correspond to this kind of states, 
if there is no symmetry breaking. So, they can be characterized by 
A 1 = ~ % aZ and A 2 = - !f z , where <r° is the 2x2 identity matrix and 
u z the usual z-Pauli matrix. The scale transformation is performed by 
A 1 = A 1 A 1 and A = A A . The Hilbert space in the block is given 
by 1 1) = 1 11) and 1 6) = 1 00). It is evident that E = E. 

3. The third case corresponds to the fixed point in the valence bond phase 
that turns out to have the same decomposition as in the fixed point of 
a one dimensional cluster state[BR00]. At the fixed point the matrices 
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appears as: A 
transfer matrix, 



iff , A x 



E= l - 
2 



Ay 



-a y A z 



^a z and the 



fl 1\ 





\1 1/ 



(4.23) 



is such that E = E 



E. 



4. The third example completes the cases where E is diagonalizable. For 
the case where E has a Jordan-block decomposition appear another 
two cases. When the local Hilbert spaces is two dimensional with 
matrices, 



A {: 







-ie > 



A 





1 



it can be seen that the A 1 is nilpotent and A°A l 
transfer matrix E is written as, 



/I 



E 



e~ l 



w 







0\ 



e ie 

1/ 



(4.24) 
' i0 A l . The 

(4.25) 



Then, the final state will be a combination of |00...010..00) vectors 
where at most one |1) appears and in the case 6 = 0, it is recovered 
the W state[DVC00]. 

5. The last case, it happens when there are domain walls. Then, the 
decomposition in terms of the matrices A s will be given by, 





cos a sin f3 e lt 
ie 



A L 




sin a 



A' 



e - 
cos a cos (3 



(4.26) 



which have a local Hilbert space of dimension three, but dimension 
two if sin a = 0. They fulfill: A A 1 = e ie A l , A 1 A = 0, A A 2 oc A 1 , 
A 2 A° 
A 2 A 2 



0, A 1 A 2 
-ie a 2 



~ ie A\ A 2 A 1 = 0, A°A° = e ie A°, A 1 A 1 



0, 



E 



e A , and the transfer matrix E will be 

/ 1 

e~ ie cos q cos P 

e cos a cos f3 

1 e td cos a sin (3 e ld cos a sin (3 I J 








(4.27) 



So, the state is a superposition of terms of the form 1 00... 00122... 22). 
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An obvious generalization of the kind of fixed point of a valence bond 
state and cluster state for D x D matrices corresponds to a symmetric su- 
perposition of D 2 orthonormal local states {\s), 1 < s < D 2 }, with the 
generators T s of J^AdxD, such that, = 1 iff £)(a — 1) + (3 = s and zero 
otherwise. Then, the local state appears as, 



D 



2 



Y^^=T s \s). (4.28) 



o=l 



with a transfer matrix 



1 D * 

e=-y,( ts )*® ts ' ( 4 - 29 ) 



which fulfills E = E 2 = E and it has just one non null eigenvector. Its 
entropy is S = log 2 D 2 independent of the number of sites. These kind of 
states are an example of fixed points in which the long range vacuum has 
non trivial entanglement properties. 



4.3 Conclusions and outlook. 

To sum up, this third part uses the matrix product state formalism to de- 
scribe any quantum state and links, at the quantum level, this formalism 
with block spin transformation that appears in statistical mechanics. We 
employ this renormalization group transformation on quantum states to 
characterize and classify the fixed points of this transformation. Then, the 
main ideas are 

1. There exists a description of real space renormalization group trans- 
formations within matrix product states. 

2. We give a complete classification of fixed points of quantum states in 
non trivial cases. 

As we have seen, this RG protocol has allowed us to obtain states that are 
fixed point of the scale transformation but have non trivial entanglement 
properties. Therefore, it is possible, within this framework, to study non 
trivial low energy sectors of a given theory [Wen04]. 



Bibliography 



[:2001] Quant. Inf. and Comp., 1, 2001. 

[Aff86] I. Affleck. Universal term in the free energy at a critical point 

and the conformal anomaly. Phys. Rev. Lett., 56:746, 1986. 

[Aff98] I. Affleck. Exact correlation amplitude for the s=l/2 heisen- 

berg antiferromagnetic chain. J. Phys. A: Math. Gen., 31:4573, 
1998. 

[AKLT87] I. Affleck, T. Kennedy, E. H. Lieb, and H. Tasaki. Rigorous re- 
sults on valence-bond ground states in antiferromagnets. Phys. 
Rev. Lett., 59:799, 1987. 

[AKLT88] I. Affleck, T. Kennedy, E. H. Lieb, and H. Tasaki. Valence bond 
ground states in isotropic quantum antiferromagnets. Com- 
mun. Math. Phys., 115:477, 1988. 

[AnsOO] D. Anselmi. Quantum irreversibility in arbitrary dimension. 
Nucl. Phys., B567:331-359, 2000. 

[Aue94] A. Auerbach. Interacting electrons and Quantum Magnetism. 
Springer, 1994. 

[BBPS96] C.H. Bennett, H.J. Bernstein, S. Popescu, and B. Schumacher. 

Concentrating partial entanglement by local operations. Phys. 
Rev., A53:2046, 1996. 

[BCN86] H.W.J. Blote, J. Cardy, and M.P. Nightingale. Conformal in- 
variance, the central charge, and universal finite-size ampli- 
tudes at criticality. Phys. Rev. Lett., 56:742, 1986. 

[BD00] C.H. Bennett and D.P. DiVincenzo. Quantum information and 

computation. Nature, 404:247, 2000. 

[Bek94] J.D. Bekenstein. Do we understand black hole entropy? 1994. 

[Ber71] V.L. Berezinskii. Sov. Phys. JETP, 32:493, 1971. 

45 



46 



Bibliography 



[Bet31] H. Bethe. On the theory of metals. 1. eigenvalues and eigen- 

functions for the linear atomic chain. Z. Phys., 71:205-226, 
1931. 

[BF64] J. C. Bonner and M. E. Fisher. Linear magnetic chains with 

anisotropic coupling. Phys. Rev., 135:A640, 1964. 

[Bha96] R. Bhatia. Matrix Analysis, volume 169. Graduate Texts in 
Mathematics. Springer- Verlag, 1996. 

[BHLM95] R.D. Ball, P.E. Haagensen, J.I. Latorre, and E. Moreno. 

Scheme independence and the exact renormalization group. 
Phys. Lett, B347:80-88, 1995. 

[BKLS86] L. Bombelli, R.K. Koul, J.H. Lee, and R.D. Sorkin. A quantum 
source of entropy for black holes. Phys. Rev., D34:373, 1986. 

[BM71] E. Barouch and B. McCoy. Statistical mechanics of the x y 
model, ii. spin-correlation functions. Phys. Rev., A3:786, 1971. 

[BPZ84] A.A. Belavin, A.M. Polyakov, and A.B. Zamolodchikov. In- 
finite conformal symmetry in two-dimensional quantum field 
theory. Nucl. Phys., pages 333-380, 1984. 

[BR00] H. J. Briegel and R. Raussendorf. Persistent entanglement in 

arrays of interacting particles. Phys. Rev. Lett., 86:910, 2000. 

[BTW02] J. Berges, N. Tetradis, and C. Wetterich. Non-perturbative 
renormalization flow in quantum field theory and statistical 
physics. Phys. Rept, 363:223-386, 2002. 

[BZJ89] E. Brezin and J. Zinn Justin, editors. Fields, Strings and Crit- 
ical Phenomena. Les Houches, Session XLIX, 1989. 

[CC04] P. Calabrese and J. Cardy. Entanglement entropy and quan- 

tum field theory. JSTAT, 0406:P002, 2004. 

[CC05] P. Calabrese and J. Cardy. Evolution of entanglement entropy 

in one-dimensional systems. 2005. 

[CDS] B. K. Chakrabarti, A. Dutta, and P. Sen. Quantum Ising 

phases and transitions in Transverse Ising models. Springer 
m.41. 

[CFL91] A. Cappelli, D. Friedan, and J.I. Latorre. C theorem and 
spectral representation. Nucl. Phys., B352:616-670, 1991. 

[CH] P. Christe and M. Henkel. Introduction to conformal invari- 

ance and its applications to critical phenomena. Spinger- Verlag 
m.16. 



Bibliography 



47 



[CH04] H. Casini and M. Huerta. A finite entanglement entropy and 

the c-theorem. Phys. Lett, B600: 142-150, 2004. 

[CHDB05] J. Calsamiglia, L. Hartmann, W. Dur, and H. J. Briegel. En- 
tanglement and decoherence in spin gases. 2005. 

[CP01] M.C. Chung and I. Peschel. Density-matrix spectra of solvable 

fermionic systems. Phys. Rev., B64:064412, 2001. 

[CW94] Jr. Callan, C.G. and F. Wilczek. On geometric entropy. Phys. 
Lett, B333:55-61, 1994. 

[DFMS97] P. Di Francesco, P. Mathieu, and D. Senechal. Conformal field 
theory. Springer, 1997. 

[dNR89] M. den Nijs and K. Rommelse. Preroughening transitions 
in crystal surfaces and valence-bond phases in quantum spin 
chains. Phys. Rev, B40:4709, 1989. 

[DVC00] W. Dur, G. Vidal, and J. I. Cirac. Three qubits can be entan- 
gled in two inequivalent ways. Phys. Rev., A62:062314, 2000. 

[ERD03] D. Esteve, J.M. Raimond, and J. Dalibard, editors. Quan- 
tum entanglement and information processing. Les Houches, 
Session LXXIX, 2003. 

[FKR04] H. Fan, V. Korepin, and V. Roychowdhury. Entanglement in 
a valence-bond-solid state. Phys. Rev. Lett, 93(22) :227203, 
2004. 

[FL98a] S. Forte and J.I. Latorre. A proof of the irreversibility of 
renormalization group flows in four dimensions. Nucl. Phys., 
B535:709-728, 1998. 

[FL98b] S. Forte and J.I. Latorre. Realization of symmetries and the 
c-theorem. 1998. 

[FNW92] M. Fannes, B. Nachtergaele, and R. F. Werner. Finitely corre- 
lated states on quantum spin chains. Commun. Math. Phys., 
144:443-490, 1992. 

[FPST94] T.M. Fiola, J. Preskill, A. Strominger, and S.P. Trivedi. Black 
hole thermodynamics and information loss in two- dimensions. 
Phys. Rev., D50:3987-4014, 1994. 

[FQS84] D. Friedan, Z. Qiu, and S.H. Shenker. Conformal invariance, 
unitarity and two-dimensional critical exponents. Phys. Rev. 
Lett, 52:1575-1578, 1984. 



48 



Bibliography 



[Fra91] E. Fradkin. Field Theories of Condensed Matter systems. 

Addison-Wesley, 1991. 

[GA89] S.M. Girvin and D.P. Arovas. Hidden topological order in 

integer quantum spin chains. Physica Scripta, T(27):156, 1989. 

[GHZ89] D.M. Greenberger, M. Home, and A. Zeilinger. Bell's The- 
orem, Quantum Theory, and Conceptions of the Universe. 
Kluwer, Dordrecht, 1989. 

[GK086] P. Goddard, A. Kent, and D.I. Olive. Unitary representations 
of the virasoro and supervirasoro algebras. Commun. Math. 
Phys., 103:105, 1986. 

[GNT98] A.O. Gogolin, A. A. Nersesyan, and A.M. Tsvelik. Bosoniza- 
tion and Strongly correlated systems. Cambridge, 1998. 

[GR94] I.S. Gradshteyn and I.M. Ryzhik. Table of integrals, series and 

products. Academic press, 1994. 

[GW02] S.D. Glazek and K.G. Wilson. Limit cycles in quantum theo- 
ries. Phys. Rev. Lett, 89:230401, 2002. 

[HLW94] C. Holzhey, F. Larsen, and F. Wilczek. Geometric and 
renormalized entropy in conformal field theory. Nucl. Phys., 
B424:443-467, 1994. 

[IJK05] A. R. Its, B.-Q. Jin, and V. E. Korepin. Entanglement in xy 
spin chain. Journal Phys. A:Math. Gen., 38:2975, 2005. 

[JK04] B. Q. Jin and V. E. Korepin. Quantum spin chain, toeplitz 

determinants and fisher-hartwig conjecture. J. Stat. Phys., 
116:79, 2004. 

[JP99] D. Jonathan and M. B. Plenio. Minimal conditions for lo- 

cal pure-state entanglement manipulation. Phys. Rev. Lett., 
83:1455, 1999. 

[Kab95] D. Kabat. Black hole entropy and entropy of entanglement. 
Nucl. Phys., B453:281-302, 1995. 

[Kat62] S. Katsura. Statistical mechanics of the anisotropic linear 
heisenberg model. Phys. Rev., 127:1508, 1962. 

[KGH+67] L. P. Kadanoff, W. Gotze, D. Hamblen, R. Hecht, E. A. S. 

Lewis, Palciauskas V. V., M. Rayl, J. Swift, D. Aspnes, and 
J. Kane. Static phenomena near critical points: Theory and 
experiment. Rev. Mod. Phys., 39:395, 1967. 



Bibliography 



49 



[KitOO] A. Kitaev. Unpaired majorana fermions in quantum wires. 

2000. 

[KM04] J. P. Keating and F. Mezzadri. Random matrix theory and 
entanglement in quantum spin chains. Commun. Math. Phys., 
252:543, 2004. 

[Kor04] V. Korepin. Universality of entropy scaling in Id gap-less mod- 
els. Phys. Rev. Lett, 92:096402, 2004. 

[KS94] D. Kabat and M. J. Strassler. A comment on entropy and 

area. Phys. Lett, B329:46-52, 1994. 

[KT73] J.M. Kosterlitz and D.J. Thouless. Ordering, metastability 

and phase transitions in two-dimensional systems. J. Phys., 
C6:1181, 1973. 

[Lau83] R. B. Laughlin. Anomalous quantum hall effect: An incom- 
pressible quantum fluid with fractionally charged excitations. 
Phys. Rev. Lett, 50:1395, 1983. 

[LLRV05] J. I. Latorre, C. A. Liitken, E. Rico, and G. Vidal. Fine- 
grained entanglement loss along renormalization group flows. 
Phys. Rev., A71:034301, 2005. 

[LMG65] H.J. Lipkin, N. Meshkov, and A.J. Glick. Nucl. Phys., 62:188, 
1965. 

[LO04] J. I. Latorre and R. Orus. Adiabatic quantum computation 

and quantum phase transitions. Phys. Rev, A69:062302, 2004. 

[LORV04] J.I. Latorre, R. Orus, E. Rico, and J. Vidal. Entanglement 
entropy in the lipkin-meshkov-glick model. 2004. 

[LR73] E. H. Lieb and M.B. Ruskai. Proof of the strong subadditivity 

of quantum-mechanical entropy. J. Math. Phys., 14:1938, 1973. 

[LRV04] J. I. Latorre, E. Rico, and G. Vidal. Ground state entangle- 
ment in quantum spin chains. Quant. Inf. and Comp., 4(1):048, 
2004. 

[LS04] A. LeClair and G. Sierra. Renormalization group limit-cycles 

and field theories for elliptic s-matrices. 2004. 

[LSM61] E.H. Lieb, T. Schultz, and D. Mattis. Two soluble models of 
an antiferromagnetic chain. Ann. Phys., 16:407-466, 1961. 

[MBF03] S. Montangero, G. Benenti, and R. Fazio. Dynamics of entan- 
glement in quantum computers with imperfections. Phys. Rev. 
Lett, 91:187901, 2003. 



50 



Bibliography 



[MN03] A. Morozov and A.J. Niemi. Can renormalization group flow 
end in a big mess? Nucl. Phys., B666:311-336, 2003. 

[M079] A.W. Marshall and I. Olkin. Inequalities: theory of majoriza- 

tion and its applications. Academic Press, 1979. 

[Mor94] T.R. Morris. Derivative expansion of the exact renormalization 
group. Phys. Lett, B329:241-248, 1994. 

[NC00] M. A. Nielsen and I. L. Chuang. Quantum computation and 

quantum information. Cambridge Univ. Press, 2000. 

[Nie99] M. A. Nielsen. Conditions for a class of entanglement trans- 

formations. Phys. Rev. Lett., 83:436, 1999. 

[NV01] M. A. Nielsen and G. Vidal. Majorization and the intercon- 

version of bipartite states. Quant. Inf. and Comp., 1:76, 2001. 

[OAFF02] A. Osterloh, L. Amico, G. Falci, and R. Fazio. Scaling of 
entanglement close to a quantum phase transition. Nature, 
416:608, 2002. 

[OL04] R. Oriis and J. I. Latorre. Universality of entanglement and 

quantum computation complexity. Phys. Rev., A69:052308, 
2004. 

[ON02] T. J. Osborne and M. A. Nielsen. Entanglement in a simple 

quantum phase transition. Phys. Rev., A66:032110, 2002. 

[OP93] M. Ohya and D. Petz. Quantum entropy and its use. Springer- 

Verlag. Berlin, 1993. 

[OR95] S. Ostlund and S. Rommer. Thermodynamic limit of density 

matrix renormalization. Phys. Rev. Lett., 75:3537, 1995. 

[Orb58] R. Orbach. Linear antiferromagnetic chain with anisotropic 
coupling. Phys. Rev., 112:309, 1958. 

[Oru05] R. Oriis. Entanglement and majorization in (l+l)-dimensional 
quantum systems. 2005. 

[PEDC05] M. B. Plenio, J. Eisert, J. Dreissig, and M. Cramer. Entropy, 
entanglement, and area: analytical results for harmonic lattice 
systems. Phys. Rev. Lett, 94:060503, 2005. 

[Pes03] I. Peschel. Calculation of reduced density matrices from cor- 

relation functions. J.Phys.A:Math.Gen., 36:L205, 2003. 

[Pes04a] I. Peschel. On the entanglement entropy for a xy spin chain. 
J. Stat. Mech., page P12005, 2004. 



Bibliography 



51 



[Pes04b] I. Peschel. On the reduced density matrix for a chain of free 
electrons. J. Stat. Mech., page P06004, 2004. 

[Pfe70] P. Pfeuty. The one-dimensional ising model with a transverse 

field. Annals of Phys., 57:79, 1970. 

[PKL99] I. Peschel, M. Kaulke, and O. Legeza. Density-matrix spectra 
for integrable models. Ann. Physik (Leipzig), 8:153, 1999. 

[PreOO] J. Preskill. Quantum information and physics: Some future 

directions. J. Mod. Opt, 47:127-137, 2000. 

[PVMDC05] M. Popp, F. Verstraete, M. A. Martin-Delgado, and J. I. Cirac. 

Localizable entanglement. Phys. Rev., A71:042306, 2005. 

[RM04] G. Refael and J. E. Moore. Entanglement entropy of random 
quantum critical points in one dimension. 2004. 

[RSW02] M.B. Ruskai, S. Szarek, and E. Werner. An analysis of 
completely-positive trace- preserving maps on m(2). Linear Al. 
and its AppL, 347:159, 2002. 

[Sch05] U. Schollwoeck. The density-matrix renormalization group. 

Rev. Mod. Phys., 77(259), 2005. 

[Sha94] R. Shankar. Renormalization-group approach to interacting 

fermions. Rev. Mod. Phys., 66:129, 1994. 

[Sre93] M. Srednicki. Entropy and area. Phys. Rev. Lett., 71:666-669, 

1993. 

[VC04] F. Verstraete and J. I. Cirac. Renormalization algorithms for 

quantum-many body systems in two and higher dimensions. 
2004. 

[VCL+05] F. Verstraete, J. I. Cirac, J. I. Latorre, E. Rico, and M. M. 

Wolf. Renormalization group transformations on quantum 
states. Phys. Rev. Lett, 94:140601, 2005. 

[VGRC04] F. Verstraete, J. J. Garcia-Ripoll, and J. I. Cirac. Matrix 
product density operators: Simulation of finite-t and dissipa- 
tive systems. Phys. Rev. Lett, 93:207204, 2004. 

[Vid99] G. Vidal. Entanglement of pure states for a single copy. Phys. 

Rev. Lett, 83:1046, 1999. 



[Vid03] G. Vidal. Efficient classical simulation of slightly entangled 

quantum computations. Phys. Rev. Lett., 91:147902, 2003. 



52 



Bibliography 



[Vid04] G. Vidal. Efficient simulation of one-dimensional quantum 

many-body systems. Phys. Rev. Lett., page 040502, 2004. 

[VLRK03] G Vidal, J.I. Latorre, E Rico, and A. Kitaev. Entanglement 
in quantum critical phenomena. Phys. Rev. Lett., 90:227902, 
2003. 

[VMDC04] F. Verstraete, M. A. Martin-Delgado, and J. I. Cirac. Di- 
verging entanglement length in gapped quantum spin systems. 
Phys. Rev. Lett, 92:087201, 2004. 

[VPC04a] F. Verstraete, M. Popp, and J. I. Cirac. Entanglement versus 
correlations in spin systems. Phys. Rev. Lett., 92:027901, 2004. 

[VPC04b] F. Verstraete, D. Porras, and J. I. Cirac. Dmrg and peri- 
odic boundary conditions: a quantum information perspective. 
Phys. Rev. Lett., 93:227205, 2004. 

A. Wehrl. General properties of entropy. Rev. Mod. Phys., 
50:221, 1978. 

X.G. Wen. Quantum order: a quantum entanglement of many 
particles. Phys. Lett, A300:175, 2002. 

X.G. Wen. Quantum Field Theory Of Many-body Systems: 
From The Origin Of Sound To An Origin Of Light And Elec- 
trons. Oxford Graduate Texts, 2004. 

S. R. White. Density matrix formulation for quantum renor- 
malization groups. Phys. Rev. Lett., 69:2863, 1992. 

K.G. Wilson. The renormalization group: Critical phenomena 
and the kondo problem. Rev. Mod. Phys., 47:773, 1975. 

K. G. Wilson and J. Kogut. The renormalization group and 
the epsilon expansion. Phys. Rep., 12:75, 1974. 

M. M. Wolf. Entropic area law for fermions. 2005. 

C.N. Yang and CP. Yang. One-dimensional chain of 
anisotropic spin spin interactions. 1. proof of bethe's hypothe- 
sis for ground state in a finite system. Phys. Rev., 150:321-327, 
1966. 



[Weh78] 
[Wen02] 
[Wen04] 

[Whi92] 

[Wil75] 

[WK74] 

[Wol05] 
[YY66] 



[Zam86] 



A. B. Zamolodchikov. 'irreversibility' of the flux of the renor- 
malization group in a 2-d field theory. JETP Lett., 43:730-732, 
1986. 



Appendix A 

Entanglement and order 
relations 



A.l Entanglement. Schmidt decomposition. Or- 
der relations 

In this appendix, we will fix our attention to the entanglement properties 
of pure quantum states in the bipartite scenario 1 . This case is particu- 
larly well understood in quantum information, see for example[:2001, NCOO]. 
In the next lines, we will follow the work of Nielsen [Nie99], mainly, and 
Vidal[Vid99, NV01] or Jonathan and Plenio[JP99], where they relate entan- 
glement properties of a quantum state to order relations. 

In a general scenario, any pure state is entangled if it cannot be written 
as the tensor product of its parts, i.e. \ip) / \4>a) ® \4>b) ® ■■■\4>n) ® 
where £ H = ®{Hi, and \4>i) G Hi- In the bipartite case, one of the 
most useful tools to describe entanglement properties is the Schmidt de- 
composition which is a consequence of the singular value decomposition for 
any complex matrix. The Schmidt decomposition says that given any pure 
state € TLa ®TLb, there is always an orthonormal basis {|<^?)} for the 
space TLa and another )} for TLb such that \ip) = J2i ) w ith 

some non-negative real numbers \ii called Schmidt coefficients that fulfill 
Si^i) 2 = 1- Then, a necessary and sufficient condition for a state to be 
entangled is that its Schmidt decomposition has more than one coefficient. 

1 A pure quantum state is described by a normalized vector, in a Hilbert space, 
TL. The bipartite scenario appears when the Hilbert space is decomposed in the tensor 
product of two different Hilbert spaces, i.e. TL = Ha <S>TLb- 



53 



54 



A. 2. Majorization relations and entropy of entanglement 



From this decomposition, a direct consequence is deduced for the reduced 
density matrix of the subsystems, which is defined by tracing out the other 
subsystem, i.e. p A = tr B \^}{^\ = Yli{Vi) 2 \4>A )(<Pa I and Pb = tr^ = 
Yli{p i i) 2 \4>B )(4>b I- ^ can be seen from their definitions that they share the 
same spectrum of eigenvalues. Then, a pure state will be separable, if the 
reduced density matrices of its parts are one dimensional projectors (pure 
states), or in other words, the state of the whole system is entangled if the 
reduced states of its subsystem are mixed. 

The amount of mixing in a density matrix can be measured by several 
mathematical tools, among them appear the entropy of entanglement and 
the majorization relations. Both quantities are deeply related and their 
connections with quantum mechanics and statistical physics are a wide and 
interdisciplinary area of research. In what follows, we will show their math- 
ematical definitions and sketch their basic properties and different meanings 
with which they can appear. 

The fact that majorization and entropy can compare the degree of mixing 
between two states imposes an order relation in the set of density matrices. 
The meaning of order is a well defined concept in mathematics: a binary 
relation -< in a set .4, in our case the set of density matrices, is called 
a preorder relation in A if and only if V x,y,z € -4: (i) x -< x and (ii) 
x -< y and y -< z then x < z. The set A with the preorder relation, i.e. 
(-4, -<) is a preordered set. A preorder relation is a partial order relation 
if (iii) given x -< y and y -< x then x = y. So, a partial order relation 
is a reflexive (condition (i)), transitive (condition (ii)) and anti-symmetric 
(condition (iii)) relation. In addition to this, a partial ordered set is totally 
ordered if V x, y G A, it happens that x -< y or y -< x. 



A. 2 Majorization relations and entropy of entan- 
glement 

Initially, majorization relations were defined in statistics and economy fields 
to compare two probability distributions [M079, Bha96]. So, given two nor- 
malized vectors x, y € lZ n , y majorizes x, written x -< y, if the following set 
of inequalities are fulfilled, 

x[ <y[ 
x[ + x l 2 <y[ + y l 2 

\ (A.l) 

x{ + ■ ■ ■ + 4_! <y{ + --- + vLi 
x[ + ■ ■ ■ + x l n =y[ + • • • + yi, 
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where the notation x\ corresponds to the components of the vector x written 
in non-increasing order, i.e. x\ > x\ Vi < j and 1 < i < n. 

As an example, it can be seen that for any normalized vector x £ lZ n , 
it always happens that (^,^,--- , ^) -< (x[,x^,--- ,xi) -< (1,0, ••• ,0). So, 
majorization coincides with the intuitive idea that when x is a more chaotic, 
disorder or mixed probability distribution than y, then x -< y. If in a set A, 
we define an equivalence relation ~ between elements x, y £ A with the same 
majorization properties, i.e. x ~ y if x -< y and y -< x then majorization 
defines a partial order in the set A. 

The connection between majorization and quantum mechanics can be 
done using the spectrum of the density matrix as a normalized probability 
distribution. In fact, given two hermitian matrices X and Y, it is defined 
X -< Y, X is majorized by Y, if A(X) -< A(Y), where A(AT) is the spectrum 
of the eigenvalues of the matrix X. 

There are other definitions for majorization equivalent to eqs.(A.l) that 
do not depend on the order of the eigenvalues of the density matrix. We 
will show three results that reveal the close relation between majorization 
and unitarity. Also, we will see that two density matrices are related by 
majorization if the more disorder matrix can be obtained from the purer 
one by a determined set of quantum operations. 

Any quantum operation £ acting on a density matrix p can be written as 
£(p) = YliEipEj. Where {Ei} is a set of operation elements that, for any 
physical evolution, they must obey the completeness relation, £V E\E{ = I 
or equivalent any quantum operation is a trace-preserving map between the 
set of the density matrices, tr(£ (p)) = tr(p) = 1. 

A specific set of quantum operation are those ones that has the identity 
matrix I as a fixed point, i.e. £ (I) = Y2i E%E^ = I. This kind of operation are 
called unital operations. The quantum operations that are unital, positive 
and trace-preserving maps are called doubly stochastic quantum operations. 
An special case of these ones are the random unitary operations, written 
£(p) = YliPiUipU} , where pi is a probability distribution and Ui are unitary 
matrices. 

The following theorem due to Uhlmann characterizes majorization in 
terms of doubly stochastic and random unitary quantum operations: 

Uhlmann's theorem - Given any pair of hermitian matrices X and Y the 
following statements are equivalent, (i) X -4 Y\ (ii) X can be obtained from 
Y applying a random unitary quantum operation £, i.e. X = S(Y); (iii) X 
can be obtained from Y applying a double stochastic quantum operation £, 
i.e. X = £{Y). 

If we write the equivalent statements for the eigenvalues of the den- 
sity matrices: (i) A(X) -< A(Y); (ii) A(X) can be obtained from A(Y) 
by a convex combination of permutations of its components, i.e. A(X) = 

pjPjA(Y), where the real numbers pj are such that < pj < 1, ^ ■ pj = 1 
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and {Pj} is a set of permutation matrices; (iii) A(X) can be obtained from 
A(Y) applying a double stochastic matrix D, i.e. A(X) = DA(Y). 

The fact that two density matrices fulfill majorization relations if they 
are related by a double stochastic operation characterizes the set of accessible 
states from a given one. This characterization is a consequence of the next 
theorem, due to Birkhoff, 

Birkhoff's theorem- The set of doubly stochastic matrices is a convex 
set whose extreme points are permutation matrices. 

So, the set {x \ x -< y, x, y G A} is the convex hull of all points obtained 
from y by permuting its coordinates. 

The third result was shown by Horn and says, 

Horn's lemma- Given two vectors x, y € lZ n with x -< y, there is an 
orthogonal matrix u such that Xi = ^ • \uij\ 2 yj. Conversely, given an unitary 
matrix u such that X{ = ^ - |it»j| 2 2/j then x ^ y. 

As we said in the introduction, apart from majorization, there are many 
other quantities that can measure the degree of mixing in a density ma- 
trix and give a notion of order in a probability distribution. Such kind of 
quantities that preserve the majorization relations are called Schur-convex 
function. A function <f> defined on lZ n is Schur-convex if given x,y £ lZ n 
such that x -< y then <f)(x) < 4>(y). 

An example of this set of function is given by any convex function, 
in particular, the entropy[OP93, Weh78] of a density matrix p defined as: 
S(p) = — trplog 2 /5. In fact, the entropy is not a Schur-convex function, but 
concave. So, given two density matrices p and a such that p^a,we know by 
Uhlmann's theorem that p = YliPiUiO~uJ '. Then S(p) = S [YliPiUiaUjj . 
Using the concave property of the entropy gives, S(p) > S(a). 

A. 3 Entanglement transformation 

Quantum information perspective has brought new insights and meanings 
to quantities like majorization or entropy. One of the new task arising for 
the studying of communication and processing of information with quantum 
mechanic support was the transformation of quantum state using local oper- 
ations and classical communication (LOCC). In a bipartite system, it is know 
that a protocol involving LOCC is equivalent to perform a measurement in 
one of the parties and applied a unitary transformation, that depends on 
the outcome of the measurement, in the second subsystem. Within this 
protocol, the question is to characterized the set of states {\4>)} that are 
accessible to a given one \ip), i.e. \tp) LOCC \(j)) . 

In 1998, Nielsen [Nie99] showed that majorization relations determined 
this set of states for a single copy of the state. His proof goes as follows; from 
the point of view of the subsystem that realizes the measurement Ej, the 
state will change to Ejp^Ej = Pjp<f>, where p^ is the reduced density matrix 
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before the measurement, p^ is the density matrix after the measurement and 
Pj is a probability constant. If we used the polar decomposition, 



E iyfH = \J E i VP^\J VN> e ] v j = y/PjP4> v j ( A - 2 ) 

with some unitary Vj. Then, p^ = J2j PjVj P<t>Vj an d 1^) 1^)- 

In 1996 Bennett, Bernstein, Popescu and Schumacher [BBPS96], using 
the singlet state as the unit of quantum correlation, showed the number m 
of singlets that can be extracted from an ensemble of n copies of a given 
bipartite state. The rate of inter-conversion is given by the entropy of entan- 
glement in the asymptotic limit. Their proof is sketched in the next lines. 
Without loss of generality, let \ip) = cos#|00) + sin(9|ll) be the Schmidt 
decomposition of a bipartite, two level quantum state, where 9 € [0, 2tt) 
is some angle and {|0), |1)} is an orthonormal local basis for every subsys- 
tem. We know that the entanglement entropy of this system is given by, 
E{\ip)) = — cos 2 9 log 2 cos 2 6 — sin 2 6 log 2 sin 2 9. If there are n copies of this 
state then, 

|V)® n = E ( cos #) n ~ N (sm0) N \x)\x) 

(A.3) 

= ^(cos#r-(sin#r E l x >l x >' 

w=0 \x\=uj 



where \x\ corresponds to the number of one's in binary representation. The 
state 

' n ' 



\ 1-1/2 

E ( A - 4 ) 

\x\=u 



is a maximally entangled state equivalent to m = log 2 ( U ) singlets. From 



uJ 



the n copies of the initial state, there is a probability p(uj) = ( cos 2 ( n -^ sin 

to obtain the maximally entangled state This probability distribution 
can be described by a gaussian distribution when n — > oo with a mean value 
nsin 2 # and variance nsin 2 #cos 2 #. So, on average and in the asymptotic 
limit, 

m = log 2 Q ~ log 2 ( ns ™ 2 - nE(\if>)), (A.5) 
where it is used the Stirling approximation, log 2 n! ~ n\og 2 n — n. 
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B.l Introduction and notation 

The behavior of two point correlation function is one of the main properties 
that characterizes a critical theory. In field theory, the propagator is defined 
by the vacuum expectation value: 

G(t,r) = (0(t,r)0(O,O)) = (0(r) e -^V(O)), (B.l) 

where in the first equality, the fields are defined in the Heisenberg picture and 
in the last one in the Schodinger picture. To study its properties, usually, it 
is performed a Wick rotation in the time direction, in which it — > r where 
t G R. So, the space-time coordinates are mapped into the complex plane 
where the real part takes into account the euclidean time direction and the 
imaginary part the coordinate direction. In these lines, we will work in this 
framework and we will note the space-time surface in complex coordinates 
as z = t + ir and z = r — ir. 

In a critical theory, two points correlations transform in a covariant way 
under a scale transformation G(r) = b~ 2x G(r /b). From this expression, it is 
assumed a covariance transformation under dilatations for a kind of operator 
called scaling operators, cp(r) = b~ x (ft(r /b), where x is the scaling dimension 
of (j). An extension to other transformations that keep the covariance of 
scaling operators (or a subset of them), shows that any local combination 
of translations, rotations and dilatations in the complex plane can still be 
considered. This set of transformations forms part of a bigger group called 
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conformal. In fact, any conformal transformation can be defined as that 
coordinate transformation which keeps the metric invariant, up to a local 
scale factor, i.e., g'^ u (r') = VL(r)g^ v (r). 

Given any infinitesimal local coordinate transformation, r'^ = r /1 +a A '(r), 
the line element defined as ds 2 = g^dr^dr" will change to ds 12 = ds 2 + 
(d^a u + d u a fjL )dr^dr v . If the transformation is conformal, the metric is 
invariant up to a constant, then {d^cty + d u a^) = d\a x gfj, u , which implies 
in the flat complex plane, where = 8n V and ds 2 = dzdz, that d z a z = 
dzOi z = 0. Therefore, the conformal group in two dimension is isomorphic to 
the group of analytic transformation in the complex plane where z — > oj(z) 
and z — ► lo{z). For example, a global translation in the plane is represented 
by z — > z' = z + a; a rotation by an angle is described by z — > z' = ze~ ld ; 
or a global dilatation can be recast as z — > z' = bz. 

Motivated by the models that appear in spin systems, as an explicit 
example, we will show how a free fermion theory looks like within this for- 
malism. The Ising model can be characterized by a hamiltonian that reads, 

with a set of anticommuting operators {St,,^} = 5^^ that evolve in the 

Heisenberg picture like, b<j> = i[H,b^\ = — or 6</,(t) = e - *'^'*^. It has 
a vacuum state |0) such that bs\0) =0 V</> and the majorana fermions 



H = I d<j> \(j>\ b\b^ = J d<P4> (&T& + ) . 



r 



J<t>>0 V 2ir v r / 



l(j>>0 

{a a ,r,a/3,r'} = S a p S(r - r'), a a , T d 



map the hamiltonian into, 

H=yJ dr (®R,rdrd R , r - a L>r d r a L>r ) . (B.4) 
The time ordered correlation, i.e. t\ > t 2 , is written as, 

(a R (h, ri )a R (t 2 ,r 2 )) = ±- f dcf> e Wn-V e -Wr*-*) ( B .5) 

in the complex coordinates where r = it with r £ R and z = r — ir, the last 
equation is recast into 

{a R { Zl ,zi)a R {z 2 ,Z2)) = [ # e-^-~^ = (B.6) 
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where 1Ze{z\} > 1Ze{z 2 }. In a similar way, it can be obtained the rest of 
time ordered correlators, 

{a R (z 1 ,z 1 )d L (z 2 ,z 2 )) = {aL(zi,zi)a R (z 2 ,z 2 )) =0, 

which yields the scaling dimensions x = \ for the operators a R and a^. Also, 
it can be obtained the correlator with the fields d z aL(z) or d z a R (z), 

{d zl d R (zi, zi)a R (z 2 , z 2 )) 

(d Zl a L ( Zl ,Zi)a L (z 2 , z 2 )) 
{d Zl a R (zi, zi)d S2 a R (z 2 ,z 2 )) 
(d zl a L (zi,Zi)d Z2 a L (z 2 ,z 2 )) 
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vr (zi - Z2 y 



B.2 Conformal symmetry 

Following Belavin, Polyakov and Zamalodchikov[BPZ84], a conformal in- 
variant theory can be defined in an axiomatical way: 

1 . The theory is described by the correlation functions of a set (in general 
infinite) of local scaling operators {4>(z, z)}. 

2. This set of operators {4>{z, z)} is supposed to be complete in the sense 
that they close the operator algebra 4>i(z)(f>j(0) = J2k Cfjiz)^^) with 
structure constants C^(z). 

3. There is a subset of operators called primary which transform under 
any conformal transformation uj(z), global or local, covariantly, 

AaM z > e )^(j£) a &aM ( b - 9 ) 

where the real numbers (A, A) are called complex scaling dimensions 
or conformal weights. 

4. Any local operator can be written as a linear combination of the pri- 
mary operators and their derivatives. 

5. The vacuum is invariant under global conformal transformations. 
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6. The generators of conformal transformations are given by the com- 
ponents of the energy- momentum tensor T^ v (z), characterized by the 
Virasoro algebra. 

Before showing in more detail these postulates, it has to be fixed the 
time direction in the complex plane. This fact is done using concentric 
circles around the origin of the complex plane as the equal time directions, 
this choice is called radial ordering. For instance, we can parametrize a finite 
model of N sites by the coordinates z = r+ir, identifying the extreme points 
in the space direction, i.e. r + N = r, so, the model is placed in a cylindrical 
geometry. Also, we can place the system in the entire complex plane using 

27T Z 2lT / \ 

the conformal transformation, uj = e~ = e^~^ T+ir ' which maps the infinite 
past r — ► — oo to the origin of the plane and the equal time directions are 
organized in concentric lines around this point. With this mapping, the time 
reversal operation r — ► — r is performed by u — > -\, and time translations 
in the cylinder r — > r + a corresponds to dilatations in the plane uj — > e a uj. 

With the postulates and the radial quantization, we will see how global 
conformal invariance fixes the functional form of the propagator between 
any two primary operator. Under any conformal transformation z — > uj(z), 
the two point function with \z\\ > \z2\ will change to 



( A A 1 ,A 1 (. Z ^ Z ^ A A 2 ,A 2 (^,Z 2 )} 



duj\ 1 f duj\ 1 f duj 



- 



dz\ ) V dz\ J \ dz2 



A 2 



(^A 1 ,A 1 (^l>^l)^A 2 ,A 2 (w2,^2)), 



(B.10) 



in the case of an infinitesimal conformal transformation, z — > z + a(z), 
primary operators transform as 

= (1 + 8 z a) A (1 + 8- z af A A>A (z + a(z),z + a(z)) - A AjA (z, z) (B.ll) 
= [(Ad z a + ad z ) + (A^a + ad z )] A A ^ A (z, z) + 0(a 2 ), 

global transformations do not change the two point function, then 

^A A ^ Ax {z x ,zx)A A ^ A2 {z2,Z2)) = 

= ( <J a>lA 1 ,A 1 ( 2 l^l)^A 2 ,A 2 (^2)} + (^A 1 ,A 1 ( 2 l) 5 «^A 2 ,A2( 2; 2^2)) 



\^ A n d z a + ad Zn + A n d z a + ad z 



(B.12) 



.71=1,2 



• (^A 1 ,A 1 (^l^l)^A 2 ,A 2 ( 2; 2,22)) = 

Translational invariance, i.e. a = 1 and a = 1, shows that the coordinate 
dependance should be like z±2 = z\ — Z2 and z\ 2 = Z\ — Z2- Scale invariance, 
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i.e. a = z and a = z, fixes the functional dependance, 

C 

(A Ai ,aM> E i) A A 2 aM>^)) = A 1+ A 2 -A 1+ A 2 ■ ( R13 ) 

Z 12 Z 12 

where C\ 2 is a constant that can be determined by normalization. Finally, 
using a = z 2 and a = z 2 implies, 

n 12 a; a" ( R14 ) 

U Ai f Zi 2 

From this expression, the scaling dimensions can be defined in terms of the 
conformal weights x = A + A and the spin of the operator 1 , s = A — A, 
which means that the two points correlator behaves as G(z, z) = \z\~ 2x (|) S . 

A conformal theory postulates that the operators A{z\) and A{z2 ) are 
related via the operator algebra, in the limit z\ — > Z2, the product of these 
operators can be characterized with the most divergent parts, so, it is usually 
defined the product expansion at the level of operators: 

C 

Z 12 Z 12 

As we saw, in the case of a free fermion Al = A# = ^ and A^ = A^ = 
and the operator product expansion (OPE) appears as, 

o.r(zi, Zi)a R (z 2 ,z 2 ) 



2tt z\ - z 2 (Tt 
1 1 (B.16) 

a L (zi,zi)a L (z2,z 2 ) ~ 

Z7T Z\ — Z 2 

Another important OPE in conformal theories is the expansion of the 
operators with the energy-momentum tensor. In the case of the fermionic 
model, the components of this tensor, in normal ordering, are written as, 

T(z) = -it : a L (z)d z a L (z) : 

= -ir lim (a L (z)d UJ a L {Lo) - {at{z)d u a L (u))) 

(B.17) 

T(z) = -7T : a R (z)d s a R (z) : 

= -7T lim {a R (z)d Q a R (uj) - {a R (z)d a R (uj))) 

z—nu 

Because this theory is gaussian, the Wick theorem can be applied in the 
product of these operators. Taking into account the contractions between 



1 Although, it is called spin, it has no relation with the spin property in quantum 
mechanic system. 
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pairs, with the sign that comes from the permutation of the fields, 

T{z)a L (uj) = -7r : a L (z)d z a L {z) : a L (oj) ~ H 

^ (z — UJ) z — UJ 

T(z)afl(a;) = -vr : 0^(2)^0^^) : or(w) ~ -— —2 + — — , 

^ [z — UJ) Z — UJ 

and the expansion of the energy momentum tensor with itself is written 

T(z)T(uj) = 7T 2 : a L (z)d z a L (z) :: a L (uj)duaL(uj) : 

1 1 27V) (B.19) 

~ 4(z-u;) 4 + (z-u;) 2 (z - u) ' 

in a similar way it can be done with the right mover or antiholomorphic 
terms. 



B.3 Virasoro algebra 

With the radial ordering in the complex plane, it is straightforward to get 
equal time commutators between two operators. For instance, given two 
local fields 4>i(z) and 02 (w), the radial ordering corresponds to 0i(z)02(w) 
if \z\ > \uj\ and 2 (u;)0i(2:) if \uj\ > z. Then, the contour integral 



dz 0i(z)0 2 (u;) ( B -20) 
is equivalent to the commutator [<f <i,20i(z), 2 (u;)], he. 



dz 0i(z)0 2 (u;) = f dz 4>i{z)<f)2{oj) - f dz 4>2(uS)4>i{z) 

\z\>\w\ J\A<\u\ /-g 21) 

cfe0i(z),0 2 (a;)] 



Figure B.l: Equivalence between equal time commutator in the complex 
plane and contour integrals. 



In the next lines, we will see that knowing the operator product expan- 
sion of the energy momentum tensor is equivalent to knowing the commu- 
tator of the operator with the generator of the conformal transformations. 

In the last section, we saw how primary operators change under a con- 
formal transformation. A way to analyze any transformation is studying its 
generators, so if an operator changes under an infinitesimal transformation, 
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this change is expressed as a commutator of the operator with the generator 
of the transformation, i.e., 5 a A A A (z, z) = [Q a , A A A (z, z)\ . For a general 
Lie group with generators gi, an infinitesimal transformation is given by e^; 
in a conformal transformation r'^ = + a^(r), because the small param- 
eter is a function of the coordinates, a general infinitesimal transformation 
is postulated as a sum in the complex plane, Q a = ^- J T^ v {r)d^ a v '{r)d 2 r , 
in analogy with the elasticity theory, which defines the energy-momentum 
tensor or stress tensor. In fact, Q a is the change in the action when an in- 
finitesimal conformal transformation is applied and obviously this expansion 
is valid only if the theory is supposed local. 

From the definition and the invariance of the vacuum under global con- 
formal transformations, it implies that: 

1. Translational invariance, i.e. a v {r) constant, is trivially obtained 

2. Rotational invariance, i.e. d^a u (r) = —d u a^(r), implies that T^ u (r) = 
T v ^(r), so the energy- momentum tensor is symmetric 

3. Scale invariance, i.e. d^a u {r) = 5^ v , made the stress tensor traceless, 

Using the fact that the energy-momentum tensor is symmetric, traceless 
and a conserved quantity, i.e. d^T^p = , it can be defined two independent 
component, 

T zz = T n -T 22 -2iT 12 , T zz =T n -T 22 + 2iT 12 , (B.22) 

where d z T zz = d z T zz = due to the energy-momentum conservations, then 
T zz = T(z) and T zz = T(z). Applying the Gauss theorem, the variation in 
the action can be recast into, 

Qa = ^~ I d"(T„Jr)a"(r))d 2 r 

= <t dza(z)T(z) + d> dza(z)T(z) 

2m J 2m J 

Finally, the variation of any primary operator can be obtained as 

= [dza(z)T(z),A AA (u,u>)] + [dza(z)f(z),A AA (u,u})] 

= 2vri / dZ a ^ z J T ^ A ^A^ u ^ u ) + 2~7ri J- dZ a ( z ) T ( z )^A,A(^ w ) 
= ((A<9 w a + adu) + (Ad Q a + acta)) A a<a (uj,u) 
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So, applying the residue theorem, 

T(z)A A {cu) ~ A . 2 A a (lj) + — 5 w A A (w) 

(z — i0) z z — UJ 

A 1 (R25) 

T(z)A a (uj,uj) ~ — — -~tA a (uj) + — — —d Q A A (uj,uj) 

\z — uoy \z — uj) 

or writing the commutation relations, 

J- [T(z),A A (u)] = 5(z - u)d u A±(u) - Ad z 5(z - co)A A (co) (B.26) 

ITTl 

and in a similar way for the antiholomorphic part. These equations match 
with the OPE we obtained for the free fermion model if Al = Ar = ^ and 
A R = A L = 0. 

A remarkable fact about the OPE of the energy-momentum tensor is 
that it does not appear as the one for primary operators, as we saw in the 
model of the free fermion. In a general case, the OPE of T(z) reads, 

T(z)T(u) ~ -7 rr + 7 —^T + , \ 

T(z)T(Q) ~ -7 rj + 7 + , r 

y ' KJ 2 (s - w) 4 (z - uf (z - u) 

where the coefficient of the first term in the expansion is a c-number that 
fixes the universality of the model, it is called the central charge or con- 
formal anomaly. In the fermion example, c = c = ~- The second term in 
the OPE tells us that the energy-momentum tensor T(z) has dimensions 
(A, A) = (2,0) and T(z) has dimensions (A, A) = (0,2). Then, the two 
point expectation value of the stress tensor is written, 

(T(z)TM) 



^"^ (B.28) 

and the variation of the stress tensor under infinitesimal conformal trans- 
formations reads, 

5 a T(to) = — I dz a(z)T{z)T{uj) 

2m J<* (B.29) 
= ^dla(uj) + 2T(u)d u a(u) + a(u)d u T(u) 

or writing it as the commutator of the energy density, 

^[T{z),T((j)]=S(z-u)d u T(u)-2d t S(z-u)T{u) + ^6(z-u) (B.30) 

AlTl 6 
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for finite transformations, it is recast into, 



f ^ = {t)' Tiz) + T2^ (R31) 

where is the Schwarzian derivative, 

(^--^sr-sl-^sr) (B - 32) 

Once we know how the energy-momemtum tensor acts on primary fields 
and their OPE, it is useful to expand the tensor and fields in Laurent series, 
such that 

= E E z- m - A z- n - A A m , n 

mdZ ndZ 

T(z) = £ z--*L„; L,=/^'T W <R33) 

n£.Z 



where the quantum generators L n obey the Virasoro algebra, 
[L n ,L m ] = (n - m)L n+m + j^n(n 2 - l)8 n+m>0 ; 

[L n , L m ] = (n - m)L n+m + ^^(n 2 - l)5 n+mfi ; (B.34) 
[L n , L m ] = 0. 

The primary fields and the L n generators are related through the following 
commutators 

= A(n + l)u n A AA (u,iD)+Lu n+1 d LlJ A AA (uj,u>) (n > -1) 

(B.35) 

= A(n + l)^M A ^(a;,w)+w ri+1 ^^ AiA (u;,a;) (n > -1) 

A direct application of these equations, and specially important in quan- 
tum spin models, is the relation between the energy density defined in the 
complex plane and the one that appears in a chain of N sites. We know 
that the transformation that maps the plane into the cylinder is given by: 

N_ ■ 
2tt 



oj = t + ir = rL log z, then 



Tcyl(u) = \z 2 T p i a {z) - ^ 



(f) E( L «-^<o) e 



2imu) /N 



(B.36) 
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which implies that (T cy i(uj)) = — ^ (t^) 2 j so the central charge can be 
measured from finite size correction in a quantum chain. Another important 
point is that the generator of the dilatations, i.e. a{z) = z, in the plane is the 
zero component in the energy momentum tensor, Lq = <f ^zT p i a (z), and we 
know, within the radial ordering, that the dilatations in the plane are related 
with time translation in the cylinder. Then, the quantum hamiltonian of a 
conformal invariant system in one dimension with N sites is written as 

M^M-^)- (B37) 

So, diagonalize the quantum hamiltonian is equivalent to find the eigenvalues 
and eigenvectors of Lq + Lq. 



B.4 Hilbert space in conformal theories 

In the last section, we have seen that the eigenstates of Lq corresponds to the 
eigenstates of the quantum hamiltonian. The vacuum in a theory is defined 
as the state with minimum energy, then, in theories where the energy is 
bounded from below, it can be fixed to zero and 

L |0> = L o |0) = 0. (B.38) 

States created by the action of a primary field into the vacuum in the infinite 
past, are called highest weight, 

|A,A> = lim A AA (r,r)|0)= lim A AA (z,z)\0), (B.39) 

T — * — OO ' 2,2— >0 

the adjoint of this state is defined by, 



(A,A| = |A,A)T = lim A AA (z,z)\0) 
= lim (0\(A AA (z,z)y = hm (0|^ AA (-,3)^^ A 

2,2— >0 ' z,z— >0 ' Z Z Z ziA z 2IX 

which gives 

(A, A|A, A) = lim J K -L (0 |A AA (i,i) J 4 AA (0,0)|0) 

z,z— >0 Z z Z Z 



(B.40) 



z 

lim uj 2A lo 2A (0\A a x(u, Q)A A A (0, 0)|0) = 1. 

U),U)— >oo 



(B.41) 



In the case of the energy tensor, 

■t 

z ,i+'* ■ - £± t — * z - 



rt (-) = E^ ?W^ = Ei^ (B.42) 
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then LL = L_ m . Prom the OPE 



00 A 1 

T(z)A A (0)= ^ z- n - 2 L n A AA (0,0)~-A A (0) + -dA A (0) (B.43) 

n=— 00 

we know that, 



Lo\A) = AIA) 



We also saw that 



[Lq, L r 



L n \A) = Vn > 0. 



(B.44) 



(B.45) 



then, if |A) is an eigenvector of Lq with eigenvalue A, L n \A) is an eigenstate 
with eigenvalue (A — n). So, as a raising operator if n < and as a 

lowering operator if n > 0. 

From every primary field, a tower of states or descendant appears, writ- 
ten 



ni,n 2 , 



;ni,n 2 , 



A, A) 



ri2 • 



.L_ ni L_ n2 ...|A,A) 



(B.46) 



and known as Verma module which are eigenstate of Lo with eigenvalue 
A + rii, and the eigenvalue of the antiholomorphic part Lq is A + fi{. 
It can happen that not all the descendant are linearly independent, then, it 
can appear null vectors. So, a final representation of the Virasoro algebra is 
given by the highest weight and their descendant removing the null vectors. 

The characterization of the null vectors will give us information about 
the unitarity in a theory. In an infinite representation, a theory is unitary if 
there is no state with negative norm. For example, we already know that, 



(A|A) = 1; (A|LiL_i|A) = 2A; (0|L 2 L_ 



a |0> = 5 ; 



(B.47) 



then, a necessary condition for unitarity is that c > and A > 0. Perform- 
ing the same analysis at the second level in the descendant, there are two 
independent vectors (L_i) 2 |A) and L_ 2 |A). Then, it appears a determinant 
that must be positive defined in a unitary theory, 



L 2 
(Li 



4A + f 
6A 



6A 

4A (2A + 1) 



which is zero iff 



A(c) 



5- C± y/(l ~ C )( 25 - C ) 

16 



(B.48) 



(B.49) 



which shows that theories with c > 25 are nonunitaries and for c < 1 it 
is possible to find theories that are nonunitary. Following in the same way, 
extracting a determinant at every level, it is possible to obtain the conditions 
that have to be fulfill a theory to be unitary. These determinants are called 
Kac determinants and they show that 
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1. If c > 1, then unitarity only constraints the values of A > 

2. If c < 1, then the set of possible unitary theories is discrete with 

c=l-V-Tv A pg (m) = ((m+ / ) f" T ^ )2 " 1 ; (B.50) 
m(m + l) MV ; 4m(m + l) v ; 

with l<q<p<m — 1 and for every value of the central charge c 
there are values of the conformal dimensions A. 

The conditions for unitariy were shown by Priedan, Qiu and Shenker[FQS84] 
and Goddard, Kent and 01ive[GK086] proved that they were necessary and 
sufficient. 



